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Chapter  1 


Everyone  who  has  seen  the  ocean  has  watched  the  waves  on  its 
surface  rise  and  fall.  A  gravitational  restoring  force  acts  on  the 
vertically  displaced  water  resulting  in  this  near-harmonic  motion.  In 
the  same  manner,  gravitational  waves  can  be  excited  within  the  ocean's 
fluid  volume  wherever  there  is  any  change  in  mass  density,  either 
continuous  or  abrupt. 

The  spatial  and  temporal  scales  associated  with  this 
"internal"  gravity  wave  motion  are  much  different  from  those 
associated  with  particle  motion  caused  by  sound.  Typical  vertical 
displacements  due  to  internal  waves  are  10-50  m  while  horizontal 
displacements  are  about  1  km  and  periods  of  oscillation  run  typically 
from  about  20  minutes  to  12  hours.1  Particle  displacements  due  to 
sound  propagation  in  the  ocean  are  dependent  on  the  source  strength, 
of  course,  but  are  generally  on  the  order  of  Angstroms  or  microns 
while  the  period  of  a  100  Hz  acoustic  signal  is  five  orders  of 
magnitude  smaller  than  the  shortest  internal  wave  periods!  The  effect 
of  the  internal  wave  motion  is  to  move  huge  quantities  of  water  very 
slowly,  thus  rearranging  somewhat  the  sound  speed  field  through  which 
acoustic  waves  must  travel. 


What,  then,  is  the  cumulative  effect  of  these  sound  speed 
fluctuations  on  the  acoustic  field?  The  quest  for  the  answer  to  this 

p 

question  really  began  in  1972  when  Garrett  and  Munk  introduced  a 
model  describing  the  statistical  properties  of  the  sound  speed  fluctu¬ 
ations  due  to  internal  waves.  This  model  was  improved  and  refined  by 
subseauent  theoretical  and  experimental  work;  for  ?  st'mn?*"/  of  *-r»ir 
work,  the  reader  is  referred  to  Vo 1 .  80  (1975)  of  the  Journal  ot 

3 

Geophysical  Research,  which  has  been  described  by  Briscoe  as  a  kind 
of  birthday  party  to  celebrate  the  field's  growing  up. 

Armed  with  a  statistical  description  of  the  medium,  a  group 
of  scientists  (the  "JASON"  group)  employed  Feynman  path-integral  tech¬ 
niques  to  describe  propagation  of  acoustic  rays  through  the  internal 

wave  field.  A  detailed  exposition  of  this  work  has  been  published  in 
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monograph  and  a  summary  has  been  presented  by  Flatte. 

The  JASON  theory  is  actually  two  theories;  one  is  valid  in 
the  region  of  very  weak  sound  speed  fluctuations  and  "short"  (~10  m) 
ranges  (the  "unsaturated"  region)  and  the  other  is  valid  in  the  region 
of  stronger  fluctuations  or  much  longer  (~10^  m)  ranges  (the  "satura¬ 
ted"  region).  Consider  acoustic  propagation  through  a  deterministic 
environment;  sound  travels  from  source  to  receiver  along  a  set  of 
eigenrays.  Now  turn  on  the  fluctuations.  In  the  unsaturated  region, 
each  of  these  equilibrium  rays  is  still  dominant  although  each  of 
these  rays  and  its  associated  ray  tube  is  perturbed  by  the  random 
sound  speed  inhomogeneities.  In  the  saturated  region,  each  ray  has 


been  split  into  many  "microrays"  created  by  the  fluctuations  and 
separated  by  more  than  a  vertical  correlation  length.  The  path- 
integral  techniques  describe  acoustic  propagation  in  each  of  these 
asymptotic  regions;  they  are  not  expected  to  apply  in  the  intermediate 
region  where  each  ray  has  begun  to  split  into  microrays  but  these 
microrays  are  still  correlated.  Some  corrections  to  the  asymptotic 
theories  have  been  made  but  the  region  of  validity  is  not  complete. 

Experiment  seems  to  validate  most  of  the  predictions  of  the 
JASON  work  in  the  regions  where  the  theory  is  expected  to  apply. 
However,  our  understanding  of  the  problem  is  still  incomplete.  As 
mentioned  above,  the  behavior  of  the  acoustic  field  at  intermediate 
ranges  is  not  totally  understood.  The  path-integral  techniques  are 
invalid  near  caustics.  Also,  anisotropy  and  inhomogeneity  in  the  sound 
speed  fluctuations  cause  the  Markov  approximation  used  in  the  theory 
to  break  down  at  some  lower  limit  on  acoustic  frequency.  Although 
measures  of  this  breakdown  frequency  (generally  around  100  Hz)  are 
available,  the  question  remains  as  to  what  significant  effect,  if  any, 
the  internal  waves  can  have  on  the  acoustic  field  at  very  low  acoustic 
frequencies. 

The  problem  of  low  frequency  (v£,  100  Hz)  acoustic 
propagation  is  particularly  important  since  measurements  of  volume 
absorption  in  the  oceanic  waveguide  typically  involve  precisely  this 
type  of  transmission  over  ranges  of  many  hundreds  of  kilometers. 

The  need  for  a  low  frequency  theory  of  acoustic  propagation 

Q 

through  an  internal  wave  field  led  some  researchers  to  con- 


sider  the  effects  of  random  sound  speed  inhomogeneities  on  acoustic 
normal  modes.  Certainly  there  are  serious  computational  disadvantages 
to  a  normal  mode  description  in  the  deep  ocean;  however,  such  a 
description  is  not  only  valid  at  low  frequencies  but  also  is  free  of 
the  problems  of  caustics  and  the  difficulties  connected  with  the 
boundary  between  the  saturated  and  unsaturated  regions  described 
above.  In  addition,  normal  modes  offer  a  framework  to  describe  energy 
transfer  into  the  ocean  bottom  while  the  ray  theory  considered  by  the 
JASON  work  does  not  allow  this  bottom  interaction.  Hence,  hoping  to 
enhance  the  results  already  obtained  through  path-integral  tech¬ 
niques,  we  turn  to  the  stochastic  differential  equations  describing 
acoustic  normal  mode  propagation  in  a  random  environment. 

The  field  of  stochastic  differential  equations  is  a 

relatively  new  realm  in  mathematics;  only  15  of  the  81  references  in 
9 

Arnold's  book.  Stochastic  Differential  Equations;  Theory  and 
Applications,  were  published  before  1960.  Therefore,  it  has  often 
been  necessary  for  physicists  to  knead  their  problems  into  somewhat 
different  forms,  with  varying  degrees  of  physical  justifiability,  in 
order  to  make  progress  using  a  recently  developed  mathematical 
arsenal.  Specifically,  the  theory  of  stochastic  partial  differential 
equations  is  not  well  developed;  true  artistry  is  often  necessary  to 
reduce  a  physical  problem  to  one  dimension. 

The  problem  of  an  acoustic  wave  propagation  through  an 


internal  wave  field  is  a  prime  example  of  an  unsavory  state  of 
affairs.  The  fact  that  we  may  consider  the  sound  speed  field  "frozen 


in  time  with  respect  to  the  acoustic  field  eliminates  the  temporal 
coordinate,  but  there  are  three  spatial  coordinates  to  contend  with. 
The  sound  speed  fluctuations  are  stationary  in  the  horizontal  plane 
but  nonstationary  in  depth  so  that  homogeneity  cannot  be  invoked  to 
reduce  the  dimensionality  of  the  problems.  Also,  the  fact  that  the 
horizontal  correlation  distance  of  the  sound  speed  fluctuations  is 
about  10  km  while  the  acoustic  wavelength  is  on  the  order  of  meters 
leads  one  to  believe  that  the  system  may  not  be  Markov! 

In  order  to  gain  some  knowledge  of  the  system,  Kohler  and 

O 

Papanicolaou  simplified  the  problem  considerably  by  investigating 
acoustic  ;.ormal  mode  propagation  of  a  single  tone  through  a 
cylindrically  symmetric  sound  speed  field.  That  is,  the  ocean  sound 
speed  knows  where  the  sound  source  has  been  deployed  in  order  to  rise 
and  fall  symmetrically  around  it.  The  sound  speed  field  was  also 
assumed  to  be  stationary  in  range.  In  this  work,  Kohler  and 

g 

Papanicolaou  laid  the  groundwork  for  the  present  study.  They 
employed  a  theory  which  they  introduced  three  years  earlier^0  in  order 
to  derive  formal  propagation  equations,  although  they  did  not  consider 
a  particular  spectrum  of  sound  speed  fluctuations  (if  they  had,  they 
might  have  found  some  trouble  arising  from  reversing  a  limit  and  an 
integral).  The  main  contribution  of  their  work  was  the  introduction  of 
rigorous  scaling  techniques  by  which  it  is  possible  to  approximate  the 
system  by  a  Markov  diffusion  process  (Appendix  C)  along  with  an 
estimate  of  the  error  based  on  the  "smallness  parameter"  e  which 
characterizes  the  strength  of  the  sound  speed  fluctuations. 


About  the  same  time,  Dozier  and  Tappert^  published  two 
landmark  papers  which  considered  the  effects  of  internal  waves  on 
time-harmonic  acoustic  normal  modes.  Using  a  nonrigorous  perturbation 
technique,  they  were  able  to  corroborate  the  formal  coupled  mode 
equations  of  Kohler  and  Papanicolaou.  In  addition,  they  found  that 
the  scattering  coefficients,  independent  of  range,  represented  con¬ 
tributions  from  the  internal  wave  spectrum  at  a  horizontal  wave  number 
equal  to  the  difference  of  two  acoustic  normal  mode  wave  numbers. 

In  the  course  of  their  work,  Dozier  and  Tappert  choose  a 
direction  of  propagation  and  integrate  the  spectrum  of  the  fluctu¬ 
ations  over  the  wave  number  transverse  to  that  direction.  They  then 
assume  that  this  modified  spectrum,  with  its  corresponding  sound  speed 
fluctuations,  can  be  used  with  a  cylindrically  symmetric  stochastic 
Helmholtz  equation  to  describe  acoustic  propagation.  Although  these 
manipulations  reduce  the  system  to  a  soluble  one-dimensional  problem, 
the  attempt  to  localize  a  normal  mode  description  is  disturbing.  In 
consequence  of  this  method  any  transverse  fluctuations  in  the  acoustic 
field  are  eliminated  not  only  on  average  but  also  for  each  realiza¬ 
tion,  even  though  the  sound  speed  fluctuations  themselves  vary 
randomly  in  azimuth. 

Crucial  to  their  derivation  of  the  coupled  mode  equations  is 
the  assumption  of  "random  phases".  That  is,  it  is  assumed  that  cor¬ 
relations  between  different  normal  modes  do  not  persist  beyond  a  cor¬ 
relation  length  of  the  sound  speed  fluctuations.  It  will  be  shown 


in  the  present  study  that  when  azimuthal  fluctuations  in  the  acoustic 
field  are  artificially  prohibited  but  random  phases  are  not  assumed, 
significant  correlation  between  different  depth  modes  can  persist  for 
many  (<  25)  correlation  lengths,  but  that  the  form  of  the  coupled 
mode  equations  obtained  by  Dozier  and  Tappert  is  nevertheless  re¬ 
covered.  In  addition,  we  use  the  full  two-dimensional  spectrum  of 
sound  speed  fluctuations  with  the  result  that  the  scattering  coef¬ 
ficients  represent  contributions  from  the  internal  wave  spectrum  at  a 
horizontal  vector  wave  number  having  a  component  in  the  radial  direc¬ 
tion  equal  to  the  difference  of  two  acoustic  normal  mode  wave  numbers. 
The  average  energy  flux  is  found  to  be  radial  at  long  ranges. 

Nevertheless,  in  spite  of  the  various  simplifications  and 
assumptions,  the  Dozier-Tappert  work  is  important  in  that  it  was  the 
first  to  consider  the  effects  of  an  internal  wave  field  on  acoustic 
normal  modes  and  to  offer  the  concept  of  resonance  scattering  in  the 
context  of  acoustic  propagation. 

A  major  consideration  of  this  dissertation  is  the  effect  of 
azimuthal  fluctuations  on  the  average  energy  flux.  For  each  realiza¬ 
tion,  the  random  depth  mode  amplitudes  are  functions  of  range  and 
azimuth.  Each  mode  amplitude  is  expressed  in  terms  of  a  set  of 
azimuthal  modes  in  order  to  reduce  the  system  to  one  dimension.  Using 
the  scaling  techniques  of  Papanicolaou  and  Kohler, ^  the  system  is 
approximated  by  a  Markov  diffusion  process  and  moments  of  the  random 
mode  amplitudes  are  calculated  using  the  corresponding  Fokker-Planck 


equation.  The  random  acoustic  pressure  is  shown  to  be  stationary  in 
the  horizontal  plane  for  an  azimuthal ly  symmetric  source.  The  coupled 
equations  for  the  autocorrelation  function  of  the  mode  amplitudes 
derived  by  suppressing  azimuthal  acoustic  fluctuations  are  still  valid 
as  long  as  each  range  function  is  interpreted  as  a  sum  over  all  the 
azimuthal  modes.  The  correlation  between  different  modes,  however, 
decreases  with  range  somewhat  more  quickly  than  in  the  azimuthal ly 
symmetric  case.  Again,  the  average  energy  flux  becomes  radial  at  long 
ranges. 

The  order  in  which  the  material  contained  here  is  presented 
is  as  follows.  Chapter  2  reviews  the  thermodynamics  of  the  ocean  and 
presents  arguments  for  including  all  of  the  internal  wave  effects  on 
acoustic  propagation  in  the  equation  of  state.  Chapter  3  discusses 
models  of  the  deterministic  and  random  components  for  the  sound  speed 
and  the  assumptions  contained  therein.  Chapter  4  reviews  determin¬ 
istic  acoustic  normal  mode  theory  and  discusses  the  acoustic  Poynting 
vector.  Chapter  5  investigates  the  effects  of  internal  waves  on  an 
(artificially)  cylindrical ly  symmetric  acoustic  field.  Many  of  the 
concepts  which  would  be  obscured  in  a  more  complex  calculation  are 
presented  here.  Chapter  6  considers  acoustic  propagation  through  a 
random  internal  wave  field  where  now  the  effects  of  azimuthal  fluctua¬ 
tions  in  the  acoustic  field  are  taken  into  account.  Chapter  7 
presents  numerical  calculation  of  scattering  coefficients  for  compar¬ 
ison  with  the  Dozier-Tappert  results.  Chapter  8  is  a  detailed  summary 
of  the  results. 


Chapter  2 


Most  discussions  of  acoustics,  and  underwater  acoustics  in 

particular,  begin  with  the  wave  equation  for  the  acoustic  pressure. 

It  is  not  clear  that  this  is  a  legitimate  procedure  in  the  present 

study  since  we  consider  two  effects,  the  acoustic  pressure  and  the 

motion  of  the  water  due  to  gravity,  or  internal,  waves,  which  may 

interact.  Of  course,  they  do  interact;  otherwise,  this  dissertation 

would  have  no  point!  However,  if  we  can  assume  that  the  scales  of 

particle  motion  caused  by  the  sound  differ  greatly  from  those  caused 

by  internal  waves,  we  may  assume  that  internal  waves  affect  sound 

propagation  only  through  the  equation  of  state.  We  therefore  begin 

our  discussion  with  the  equations  of  motion  for  ocean  water  particles 

due  to  the  combined  sound  and  internal  wave  fields.  Much  of  the 

11  12 

following  thermodynamic  description  follows  Eckart.  ’ 

Seawater  is  here  treated  as  a  solution  of  a  single  compound 
in  pure  water.  The  internal  energy  z  as  a  function  of  specific 
volume  v,  entropy  n  ,  and  salinity  S  can  be  used  to  calculate  the 
pressure  p,  absolute  temperature  T,  and  chemical  potential  u;  i.e., 


Small  increments  of  these  quantities  are,  then. 


Some  elements  of  the  square  matrix  in  Eq.  (2.2)  are  related 
to  familiar  response  functions: 


where 


Y  =  £  (y-D 
Z  =  T/CyS 


(2.3) 


p  =  1 / v  =  mass  density. 


2  /  5  p  \ 

c  =  l-r-)  defines  the  sound  speed, 
v  3p  'nS 

CvS  =  specific  heat  at  constant  volume  and  salinity, 
Cps  =  specific  heat  at  constant  pressure  and  salinity, 

y =  Cps/Cvs  *  and 


1  (It) 


a 


=  coefficient  of  thermal  exoansion. 


D/Dt  =  3/^t+u*V  =  total  time,  or  convective,  derivative, 
u.  =  particle  velocity, 
g  =  gk,  the  gravitational  acceleration, 
f  =  twice  the  earth's  angular  velocity  of  rotation, 

A  =  stress  tensor  of  molecular  viscosity  (for  an  explicit 
definition  see  Refs.  12  and  13), 

£  =  result  of  all  other  forces  such  as,  for  example,  an 
oscillating  piston, 

G  =  the  heat  generated  by  irreversible  processes, 
h  =  heat  flux,  and 
s^  =  salinity  flux. 

Note  that  k  points  vertically  downward. 

If  we  assume  that  the  particle  motion  is  small,  we  may 
linearize  these  equations  by  writing  each  quantity  as  a  sum  of  its 
time-dependent  and  its  static  values.  Static  values  will  be 
subscripted  with  zeroes;  time-dependent  quantities  will  be  subscripted 
with  ones.  Therefore,  P=P0+Pi>  u.=u^,  f.=£.i»  p=po+cV  etc* 

Equations  (2.6)  are  ranked  in  successive  approximations,  the  zero- 
order  set  being  the  case  when  no  time-dependent  quantities  are 
considered,  the  first-order  set  where  products  of  zero-  and  first- 
order  quantities  are  considered  but  not  products  of  first-order 
quantities,  and  so  on.  Relations  between  zero-order  quantities  may 
be  exploited  in  the  first-order  equations,  thus  preserving  the  ranking 
of  the  approximation. 


In  order  to  obtain  the  zero-order  relations,  we  neglect  A 
and  all  time  derivatives.  We  assume  that  S  does  not  vary  with 
latitude;  thus,  all  zero-order  quantities  depend  only  on  depth.  Thus 


vp0  =  p0a 


and 


!!o  ,  c2  !!o  .  v 
3 z  o  3z  ro  3z 


+  K 


O  3z 


Combining  these,  we  define 


(2.7) 


(2.8) 


The  quantity  N(z),  called  the  Brunt-Vaisaia  or  buoyancy 

frequency,  has  an  important  physical  meaning.  Consider  a  small  mass 

of  fluid  displaced  adiabatical ly  a  small  vertical  distance  from  its 

zero-order  position  and  then  allowed  to  move  freely.  The  fluid 

element  will  oscillate  about  its  original  position  with  frequency  N, 

? 

provided  that  N  is  real.  If  N  is  negative,  the  fluid  motion  is 
unstable;  a  slight  perturbation  will  cause  wide  departures  from  the 
fluid  element's  original  position. 14-16 

It  is  very  useful  to  consider  the  concept  of  a  "potential 
gradient";  that  is,  the  measured  gradient  minus  what  the  gradient 
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would  be  in  a  static,  isohaline,  isentropic  ocean.  These  gradients 


will  be  subscripted  with  P.  For  example. 


j_/V i  j_!!°  a 

P0'32/p  po  92  c2 


(2.10) 


Thus,  it  is  the  potential  density  gradient  which  determines  the 
ocean's  static  stability,  i.e., 


Similarly 


(2.11) 


/  9T  \  9T  /  9T  \  9T 

(t?)  ■ -sM-sr)  .  •  (2-12) 


n  S 
o  o 


a  c 
o  o 


The  first-order  linearized  hydrodynamic  equations  are 


1 


9u, 

po  +  Vp1  +  P0  jf xu_-j  =  +  f.  +  »  (2.13a) 


3P1 

TtT 


+  v  •  ( oQ  Uj, )  =  o 


(2.13b) 


3n-i  1 


(2.13c) 


9S, 

ir  + v7S0 


—  7*  S 
°0  ' 


(2.13d) 
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with  appropriate  boundary  conditions.  In  many  cases,  it  is  more 
convenient  to  employ  the  following  form  of  Eqs.  (2.13): 


where 


3U 

3=  +  c0  [Wlc]p  +  N2(z)i<W  +  fxU  =  cf 

H  +  co  O'1'*]'  -  =  "t 


U  =  H,<d0c05,/2 
n2U)w  -  (p0c0)'/2 


/C  0/2 

£-k)  (V'=-)  ' 

r  =  r  (poco)‘1  97  (ooco>  -  “V 

c 


» 


(o0c 


o>"1/2 (irr)  ,  ■ 


Boundary  conditions  are  now 


n-U  =  0  (2.15a) 

on  rigid  surfaces  and 

=  (2.15b) 

at  the  surface  of  the  ocean.  We  see  that  when  ^  =0,  W(pqcq  )"^  is 
simply  the  vertical  displacement  of  the  fluid  particle. 

If  we  consider  a  region  of  the  ocean  far  from  heat  and 
pressure  sources,  we  may  take  ^  and  in  Eqs.  (2.14-2.15)  to  be 
zero.  Then,  the  solutions  of  these  equations  may  be  time  harmonic; 
time  derivatives  are  replaced  by  the  factor  -iw,  where  cj  is  the 
angular  frequency.  Actually,  a  general  limited  solution  to  these 
equations  as  they  stand  will  be  a  superposition  of  these  time  harmonic 
solutions  of  the  free  equations. 

It  should  be  remembered  that  P,  W,  and  U_  are  first-order 
quantities;  r  and  N  are  zero-order  quantities.  Therefore,  to  first 
order,  it  seems  that  the  hydrodynamic  equations  should  decouple  into  a 
set  describing  only  the  internal  waves  and  another  set  describing  the 
acoustic  propagation. 

To  support  this  assertion,  we  consider  an  instructively 
simple  situation  where  S,  c  ,  and  N  are  constant  everywhere  and  r  is 
zero.  This  situation  is  physically  impossible  for  the  ocean  but  not 


I 

so  bad  for  the  atmosphere  (which,  of  course,  doesn't  have  salinity). 
Our  source  terms  are  also  taken  to  be  zero.  We  do  not  neglect 
|  rotation,  but  we  do  specialize  it  by  making  our  system  a  flat  layer  of 

thickness  H  rotating  about  a  vertical  axis  with  angular  velocity 
f=-fi<.  The  solutions  to  the  field  Eqs.  (2.14)  are 

* 

P  =  FP] (z)  e‘iwt  ,  (2.16a) 

-  W  =  FW-|  (z)  e”ia)t  ,  (2.16b) 

where 

m 

i 


-y  +  —  a-  +  a  [F  ~  0 
,3xz  dyc 


(2.17) 


(a  is  the  separation  constant)  and 


dP-i  (z)  o  o 
—35—=  (N  -</ ^(z) 


(2.18a) 


dW-j  (z) 
'0  dz 


^2c2 

1  ■  ypi(z) 

03  “  T 


(2.18b) 


with  boundary  conditions  -PcQ=gW  at  z=0,  and  W=0  at  z=H.  Then, 


W1  =  w0  s^n  [f?<z-H)] 


(2.19a) 


P1  =  po  C0S[^(Z-H)] 


(2.19b) 


0'.Y.y.Y, 


ifcY.'i.Y.Y.' 


NMiMMii 
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where  6,  the  vertical  wave  number,  is  related  to  frequency  co  and  to 
horizontal  wave  number  a  by 

g  Cq(oj  -f  )  +  (N  -w  ) (oj  -f  -a  CQ)  -  0  .  (2.20; 

On  the  other  hand,  the  free  surface  boundary  condition  gives 

2  2 

N  -oi  =  g6  tanBH  .  (2.2V 

A  graph  of  the  simultaneous  solution  to  Eqs.  (2.20)  and 

17  2  2  2  2 

(2.21)  is  shown  in  Fig.  1.  Case  A  is  N  <a  c  +f‘ ;  case  B  is 

2  2  2  2 

N  >a  cQ+f  .  in  either  case,  each  solution  to  Eq.  (2.21)  intersects 
the  solution  to  Eq.  (2.20)  twice,  resulting  in  two  sets  of  normal 
modes  for  the  system.  The  dispersion  relation  for  these  modes  is 
shown  in  Fig.  2.*7  The  gravitational  (internal  wave)  modes  are 
confined  between  the  rotational  and  buoyancy  frequencies  whereas  the 
acoustic  modes  have  frequencies  much  higher  than  N.  For  the  sake  of 
clarity,  the  spacing  of  the  acoustic  modes  has  been  greatly  reduced 
and  the  slopes  of  the  other  curves  greatly  increased  in  Fig.  2. 

What  happens  in  the  ocean  is  this:  internal  waves  move  great 
quantities  of  water  up  and  down  (and  also  back  and  forth,  but  this  is 
irrelevant  to  sound  propagation)  at  velocities  slow  enough  to  essen¬ 
tially  maintain  the  static  entropy  and  salinity  gradients.  The  phase 
velocities  of  these  internal  waves  are  sufficiently  smaller  than  the 
sound  speed  that  this  motion  can  be  considered  incompressible;  the 


frequencies  are  so  low  as  to  make  gravitation  and  the  earth's  rotation 
important.  Sound  waves  are  transmitted  along  at  velocities  so  high 
that  particles  under  the  influence  of  internal  waves  seem  "frozen". 

The  sound  frequencies  are  so  high  that  gravitation  and  Coriolis 
effects  affect  sound  only  in  the  fact  that  the  medium  through  which  it 
is  propagating  has  been  rearranged  somewhat.  Also,  the  heat  and  salt 
contained  in  a  volume  of  water  displaced  by  a  sound  wave  don't  have 
time  to  diffuse  into  the  new  environment  before  the  volume  element  is 
pushed  back  again.  In  the  linear  regime,  sound  is  assumed  to  consist 
strictly  of  compression  waves. 

The  total  particle  motion  is  therefore  the  sum  of  rota¬ 
tional,  incompressible  motion  due  to  the  internal  waves  and 
irrotational ,  compressible  motion  due  to  sound.  Each  of  these  effects 
is  assumed  to  independently  satisfy  the  linearized  hydrodynamic  equa¬ 
tions;  the  effect  of  the  internal  waves  upon  the  sound  is  manifested 
through  the  equation  of  state.  That  is,  the  sound  speed  is  considered 
to  possess  a  small  perturbation  due  to  the  displacement  of  isodensity 
surfaces  by  the  internal  waves. 

We  may  now  feel  comfortable  in  deriving  the  acoustic  wave 
equation  in  the  usual  manner,  secure  in  believing  internal  wave 
effects  to  be  contained  in  the  sound  speed,  with  the  result 

2 

V^p  -  \  — w  =  source  terms.  (2.22) 

at 


Chapter  3 


The  ocean  sound  speed  is  a  function  of  horizontal  position 
r,  depth  z,  and  time  t.  We  take  as  our  model 

c(r,z,t)  =  c(z)  +  5c(r,z,t)  ,  (3.1) 

where  F(z)  is  the  mean  sound  speed  profile  and  5c(£,z,t)  describes  the 
sound  speed  fluctuations  due  to  internal  waves.  In  general,  c(z)  will 
vary  with  horizontal  position  and  time.  However,  these  changes  are 
usually  very  slow  compared  with  the  rate  of  variation  in  6c;  c(z) 
mainly  depends  on  depth. 

A.  The  Mean  (Deterministic)  Sound  Speed  Profile 

The  mean  sound  speed  profile  for  the  deep  ocean  is 
characterized  by  a  sound  channel.  Oue  to  a  negative  temperature 
gradient  c"(z)  decreases  as  depth  from  the  surface  increases.  This 
effect  is  eventually  balanced  by  the  weight  of  the  water  and  c(z) 
reaches  a  minimum  cQ  at  the  channel  axis  depth  zQ.  For  z>zQ,  the 
sound  speed  increases  with  depth  due  to  increasing  pressure  in  the 
nearly  isothermal  deep  ocean. 

Details  of  this  picture  vary  from  ocean  to  ocean  (e.g., 
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Urick  ).  For  example  cq  may  be  found  very  near  the  surface  in  cold 
Arctic  regions,  resulting  in  a  sound  speed  profile  which  essentially 
increases  throughout  the  entire  water  layer.  Nearer  the  equator,  in 
cloudy,  windy  environments,  turbulent  mixing  may  result  in  an 


isothermal  layer  just  below  the  ocean  surface.  Within  this  layer,  the 

sound  speed  gradient  vanishes  twice,  once  near  the  surface  at  the  base 

of  the  "mixed  layer"  and  again  deep  in -the  ocean  at  depth  z  . 

We  need  not  consider  these  details  for  the  study  at  hand. 

We  do  require  a  model  of  "c(z)  which  reasonably  describes  the  sound 

speed  profile  for  much  of  the  world's  deep  ocean  environment  and  has 

its  basis  in  physical  reality.  Therefore,  the  "canonical"  sound  speed 
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profile  derived  by  Munk  is  chosen. 

Munk's  derivation  of  the  sound  speed  profile  is  reproduced 

in  Appendix  A.  Briefly,  realistic  assumed  dependences  on  temperature 

and  salinity  by  the  sound  speed  are  combined  with  an  exponentially 

decreasing  stratification  profile,  which  is  "perhaps  the  most 
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intrinsic  property  of  the  abyssal  oceans"  (Fig.  3,  Munk;  also. 

Refs.  2,  20-21).  By  "stratification"  is  meant  only  that  the  density 
changes  with  depth.  The  word  does  not  always  imply  layers  but  applies 
to  continuous  density  gradients  as  well.  Although  the  assumption  of 

exponential  stratification  is  really  invalid  near  the  surface  (see 
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references  cited  above;  also  Eckart  )  we  follow  Munk  and  everybody 
else  in  ignoring  this.  The  fluid  is  statically  stable  when  the 
buoyancy  frequency 


N(z)  2 


-z/B 


(3.2) 


is  real.  In  Eq,  (3.2),  B  is  the  scale 
is  the  acceleration  due  to  gravity,  and 
density  gradient. 


depth  of  the  stratification, 
the  potential 


g 
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N  in  cycles/hour 


FIGURE  3 

THE  VAISALA  (OR  BUOYANCY)  FREQUENCY  N(z)  =  [gp'1  9 zpp]  ’/2. 
HORIZONTAL  BARS  GIVE  RANGE  IN  VALUES  MEASURED  OVER  THE  BISCAY 
ABYSSAL  PLAINS  (46°N  8°W)  BY  PINGREE  AND  MORRISON.21  BENEATH  THE 
THERMOCLINES  THESE  HAVE  BEEN  FITTED  BY  N(z)  =  NQe  Z/B,  WITH  NQ  =  2.8  cph 
DESIGNATING  THE  "SURFACE-EXTRAPOLATED"  VALUE  OF  N(z),  AND 
B  =  1.65  km  A  "BUOYANCY  SCALE."  IN  THE  PACIFIC  OCEAN  (CIRCLES)  THE 
THERMOCLINE  IS  MUCH  SHALLOWER,  TYPICALLY  100  m,  AND  A  "SUPER 
EXPOTENTIAL"  RISE  TAKES  PLACE  BETWEEN  800  m  AND  THE  THERMOCLINES. 

(FIGURE  TAKEN  FROM  REF.  19). 
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The  sound  speed  channel  appears  quite  naturally  at  a  depth 
very  near  the  scale  depth  of  the  buoyancy  frequency.  The  analytical 
expression  for  <f(z)  is 


c(z)  =  c(zQ)  exp 


(3.3) 


I 


where 


(3.4) 


and 

e  =  \  bYa 

* 

I 


(3.5) 


The  quantity  y A  is  the  fractional  sound  velocity  gradient  in 
an  adiabatic  isohaline  ocean.  For  real  oceans. 


z)  ~  c 


(3.6) 


which  is  the  well-known  Munk  (canonical)  sound  speed  profile. 

B.  The  Random  Sound  Speed  Perturbation 

The  sound  speed  at  position  x  (as  measured  in  a  three- 


dimensional  coordinate  system)  and  time  t  is  assumed  to  vary  from  its 
mean  value  "c(z)  by  the  slow  advection  of  an  isodensity  surface  by 
internal  waves.  The  advection  is  slow  enough  to  be  quasistatic; 


therefore  we  may  write  the  sound  speed  fluctuation  at  position  x  and 


where  q  is  the  vertical  particle  displacement  due  to  internal  waves. 
Our  problem  is  now  to  describe;;.  There  are  two  popular  ways  of  doing 
this.  One  method  is  to  make  whatever  assumptions  are  necessary  to 
make  the  hydrodynamic  equations  describing  the  internal  wave  field 
tractable  and  use  the  resulting  theoretical  description.  The  second 
is  to  take  the  description  of  the  internal  wave  field  directly  from 
the  experimental  results  without  explaining  in  detail  the  mechanisms 
causing  the  internal  wave  motion.  In  the  end,  we  shall  choose  the 
second  approach,  but  it  is  worthwhile  to  consider  the  first.  The 
usual  approach  is  to  assume  that  the  hydrodynamic  equations  may  be 
linearized  and  to  neglect  shear  coupling.  The  experimental  evidence 
concerning  these  assumptions  will  be  discussed  later.  The  linearized 
equations  of  motion  are 


where  the  subscript  i  denotes  internal  wave  (as  opposed  to  acoustic) 


effects  and  where  w.=£  is  the  vertical  component  of  u ■ ,  the  particle 

velocity  due  to  the  internal  waves.  Given  Eqs.  (3.8),  one  can  show, 

2  2? 
for  ir»g  /c*. 


(3.9) 

2  /  3^  . 

where  V„  =  [ — +  — *-)  (see  Appendix  B).  At  this  point  we  must  make 
n  \3x^  3y  / 

some  approximations  in  order  to  make  this  equation  tractable.  The 
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"traditional  approximation",  so  called  by  Eckart  because  all 
previous  treatments  of  the  problem  agreed  to  make  it,  consists  of 
neglecting  products  of  w^  and  f^,  the  horizontal  component  of  f.  This 
approximation  is  often  justifiable.  However,  as  Eckart  points  out, 
neglect  of  this  term  without  also  neglecting  similar  products  of  f^ 
and  the  other  velocity  components  robs  the  equations  of  motion  of 
their  self-adjointness.  The  nature  of  the  field  equations  is  thus 
completely  altered  and  radically  different  phenomena  will  be 
described.  Therefore,  if  w.f^  is  to  be  neglected,  so  must  the  other 
terms  involving  f^.  This  is  clearly  invalid  near  the  equator. 

However,  if  we  adopt  the  traditional  approximation, 


(V2wi )  +  N2(z)(v2w^ 


r2  3% 


-N2(z)  32  f2 

g  ^7  3z  rz  3z 


Plane  wave  solutions  to  Eg.  (3.10)  have  the  form 


v  «(z) 


where  (dropping  the  subscript  on  f) 


(3.10) 


(3.11) 


£mui+„2  iw-r  i.  H(t)  .  <3., 2) 

dz‘  -  T  g  32 


The  right  side  of  Eq.  (3.12)  may  be  neglected  since  the  ratio  of  the 

first  to  the  second  derivative  of  W(z)  is  no  larger  than  the  depth  of 
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the  ocean  and 


H  <  5*10‘4 


(3.13) 


The  quantity  W(z)  is  subject  to  boundary  conditions  at  the 
surface  (z=0)  and  at  the  ocean  floor  (z=H).^  W(H)  is  generally 


•  •  o  o  .  •  - .  •  .  •  .  • 
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taken  to  be  zero  since  the  ocean  floor  is  considered  to  be  imperme¬ 
able.  Also,  since  the  internal  waves  in  the  ocean  produce  only  very 

pc 

small  vertical  displacements  of  the  free  surface,  W(0)  is  generally 
taken  to  be  zero,  too. 

Given  the  boundary  conditions  and  a  particular  horizontal 

wave  number  a,  the  solutions  of  Eq.  (3.10)  form  a  discrete  set  of 

eigenfunctions,  each  corresponding  to  an  eigenfrequency  Cl  .  .  For  real 

Ja 

a,  these  eigenfrequencies  will  lie  in  the  interval  f<Jija<  max  N(z). 
A  vertically  displaced  fluid  element  will  oscillate  vertically  with 
the  buoyancy  frequency  N(z).  If  the  displacement  is  not  vertical,  the 
restoring  force  is  less  and  the  fluid  element's  motion  will  be 
elliptical  with  a  reduced  frequency.  At  frequencies  just  slightly 
higher  than  f,  the  particles  will  travel  in  nearly  horizontal, 
circular  orbits,  the  velocity  vector  rotating  anticyclonical ly 

i  pc 

(Fig.  4;  for  more  detail  see  Munk  and  Phillips  ). 

For  an  exponentially  decreasing  buoyancy  frequency 
N(z)=N  e‘z^B,  the  solutions  to  Eq.  (3.10)  are  exact,2  and  given  by 


INERTIAL  FREQUENCY 


BUOYANT  FREQUENCY 


FIGURE  4 

THE  WAVENUMBER  VECTOR  K  =  (a,  t3),  GROUP  VELOCITY  cg  AND  THE  HODOGRAPH 
OF  THE  PARTICLE  VELOCITY  u,(t)  NEAR  INERTIAL  FREQUENCY  (ft  =  f+)  AND 
BUOYANT  FREQUENCY  (ft  =  N-),  RESPECTIVELY.  K  IS  NORMAL  TO  BOTH  cg  AND  u, 

(FIGURE  TAKEN  FROM  REF  1) 
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is  a  normalization  constant  and  J  ,Y  are  Bessel  functions  of  the 

3  a 

first  and  second  kind.  The  eigenfrequencies  are  determined  by 

JOt 

the  upper  boundary  condition  Wj(a,z)=0  at  z=0. 

Now  we  may  represent  the  vertical  displacement  of  an 
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isodensity  surface  by  a  superposition  of  internal  wave  modes, 

V*  f  i(a*r-Q-  t) 

C(x.t)  =2  daGj (a)Wj (a,z)  e  Ja  ,  (3.17) 

j  J 

where  the  quantities  Gj(a)  are  complex  Gaussian  random  variables 
satisfying 

<  G  .  (a) >  =  0 

J 

<Gj(a.)G*,(a')>  =  P.  (a)6(a-a'  )6  j  ■ ,  .  (3.18) 

Before  discussing  the  form  of  the  spectrum  Pj(a)>  it  is 
necessary  to  say  a  few  words  about  the  above  analysis.  The  derivation 
of  Eq.  (3.9)  neglects  several  non-neg 1 igible  effects:  forcing,  shear 
dissipation,  nonlinear  coupling.28-29  It  is  doubtful  that  well- 
defined  internal  wave  modes,  other  than  perhaps  the  gravest  mode,  will 
have  time  to  establish  themselves  since  vertical  propagation  times  of 
the  wave  field  are  comparable  to  typical  interaction  times.  This 
suggests  that  it  might  be  better  to  replace  the  discrete  sum  in 
Eq.  (3.17)  by  an  integral  over  a  continuous  frequency  range. 


Observations 


tend  to  support  this  latter  notion;  propagating 
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waves  dominate  at  low  frequencies,  standing  modes  may  dominate  at  the 
highest  frequencies  (j~l)»and  both  descriptions  are  adequate  in 
between.  We  therefore  choose  the  second  method  to  describe  the 
internal  wave  field,  specifically  the  statistical  description  of  C, 
relying  heavily  on  experimental  work  which  will  be  referenced  as  we  go 
along.  In  this  latter  description,  one  assumes  a  random  state  with  a 
continuum  of  energy  in  frequency/wave  number  space.  Here,  a 
difficulty  arises  in  describing  the  spectrum  since  the  medium  is  not 
depth  homogeneous.  Nevertheless,  the  quantity 


can  be  considered  a  local  vertical  wave  number,  in  the  WKBJ  sense. 

Thus,  for  n<N(z),  we  can  use  Eq.  (3.15)  as  the  frequency/wave  number 
dispersion  relation. 

The  continuum  description  used  here  was  introduced  by  Garrett 
and  Munk,^’"^  modified  by  experiment, 32*34,35  ancj  stream] -jnec|  by  the 

Or: 

analytical  work  of  Desaubies. 

Any  function  of  two  (three-dimensional)  positions  x=(r,z) 
and  jc'  =  ( r^'  ,z)  relative  to  some  fixed  coordinate  system  can  be  written 
in  terms  of  their  average  position  (x'+x)/2  and  their  difference  in 
position  A£=(x-x'  ) .  The  covariance  <  r,  (r  ,z)r,  (r 1  ,z ' )  >  will  depend  on 
the  average  position  only  through  the  average  vertical  position 


n=(z+z')/2  (Ref.  37).  Thus,  if  (a,B)  is  the  wave  number  vector 
conjugate  to  (r-r_*  ,  ^),  where  ^=(z-z'(»  we  may  write 

'i— ’ (— ”— ' )  ,  (3.20) 


;(r,zk(r'  ,z‘)>  =  f  d3  J" daFr(a,?;n) cos  e1 


where  F;.(a,B;n)  is  normalized  such  that 


J  dad6F^(ot.,S;n)  =  <?2(x-x' ,n)>  =  <£2(n)>  .  (3.21) 


If  F^(a,6;n)  depends  on  the  horizontal  wave  number  ^  only 
through  its  magnitude  (that  is,  if  the  medium  is  statistically 
isotropic  in  the  horizontal  plane),  we  introduce 


P  (a»8;n)  =  2iraF;,(a,B;n) 


(3.22) 


such  that 


C  °°  r  oo  r2-rr  r  »  r 

J  d SJ  da  P.(a,3;n)  =  J  d dj  d3  I 


daaFr(a, 3;n)  -  (n)>  . 


(3.23) 

If  x=(r,<p,z)  d  _x' =(r  1  ,cp*  ,z 1 )  are  two  positions  in  some 
fixed  (cylindrical)  coordinate  system,  we  may  write  Eq.  (3.20)  in  what 
will  prove  to  be  a  more  convenient  form: 


<r(r,zk(r' , z ' ) > 


d2J  da|2  dnP^(l’°; 


;rO  cos  C- 


f 


x  e 


ia[rcos(^-^)-r'cos(<?l-9)  ] 
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(3.24) 


Finally,  we  may  transform  from  (a, 8)  space  to  (a,,Q)  space 
by  use  of  the  dispersion  relation  (3.19): 


P^(a,.Q;n)  =  Pr (a,8(a,ft) ;n)  |j|  =  MA(a)0(&) 

where,  according  to  the  references  cited  above  Eq.  (3.20), 


(3.25) 


r 

jk 


A(a)  =  -p  2 
a  +  a* 


(3.25a) 


0(n)  •  ^ 

n3 


(3.25b) 


and 


M  =  —  <c2(n)> 


(3.25c) 


ensure  that 


fN(n)  f 

J  daJ 


dn  |  daP(a,Q;n)  =  <  C2(n)> 

o 


(3.26) 


to  order  f /N(p.) •  In  most  cases  of  interest,  min  N(n)»f.  The  quantity 
t*=t(C2-f2)1/2  is  a  number  character istic  of  the  bandwidth;  t  is  an 
experimentally  determined  constant  on  the  order  of  10  -  10  h/m. 

The  model  presented  in  Eqs.  (3.25)  depends  essentially  on  a 
WKBJ  approximate  solution  to  Eq.  (3.12)  and  is  therefore  not  very  good 
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for  frequencies  near  N(n)  (Fig.  5;  taken  from  Ref.  35).  In 
particular,  the  spectrum  shows  a  pronounced  hump  just  below  the 
buoyancy  frequency,  followed  by  a  sharp  cutoff. 

This  hump  is  also  found  in  the  internal  wave  experiment 

(IWEX)  data  ’  and  in  previous  measurements  by  Voorhis  and 

Gould.'1  Vertically  propagating  waves  with  frequency  are  reflected 

at  the  turning  depth  (where  2-  N(z));  each  of  these  wavefunctions  has 

therefore  an  inflection  point  at  that  depth.  Waves  which 

destructively  interfere  in  shallower  water  are  therefore  locked  in 

40 

phase  near  the  turning  depth.  A  Langer  method  of  treating  the 
internal  wave  field40, 4^  has  been  shown  to  describe  the  observed 
frequency  spectrum  very  well  (Fig.  6).  However,  most  of  the  energy 
in  the  internal  waves  is  contained  in  frequencies  near  f  (Eq.  3.24b 
has  a  sharp  peak  at  s2=f/T75).  Thus  the  simpler  WKBJ  description  is 
adequate  for  our  purposes. 

The  description  is  not  yet  complete;  we  still  have  to 
2  2 

specify  <  ^  (n)>.  Garrett  and  Munk  took  the  WKBJ  solution  to 
Eq.  (3.14)  and  "depth  averaged  over  many  wiggles"  to  find  that,  to 
some  approximation,  <  '  (n)><*N  (n).  Whether  or  not  this  was  a 

legitimate  procedure,  the  observations  support  this  result.  When  the 
autospectra  were  measured  at  different  depths  (data  taken  in  IWEX)  and 
the  "WKBJ  normalized",  i.e.,  multiplied  by  N(z)/N(zr)  where  is  a 
reference  depth,  the  variation  in  the  spectral  levels  were  reduced  to 


ENERGY  DENSITY 
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some  10%32’34  (Fig.  7)  except  near  the  buoyancy  frequency,  where  the 
WKBJ  approximation  breaks  down.  Again,  we  ignore  the  region  of 
invalidity. 

Our  model  of  the  sound  speed  fluctuations  is  now  finished. 
Relating  these  sound  speed  fluctuations  to  the  displacement 
correlations  gives 


<<Sc(r,<?,z)6c(r'  ,(p'  ,z* )  > 
c(z)c(z' ) 


N  f  N2(z)  N2(z‘ ) 
V  NTnl 


f  N(n)  r 

Jf  dQl 


,  a  ,_2  f2J/2 

dfi  I  da  — g-'V"  - T - cos(ct^X(p)) 

(a  +a*)  PJ 


-  ^  j  ^ 

|  de[e1a^rcos^_0)*r,cos(<P'-e)]  +  eia[rcos(^+0)-r'cos(<A'  +  e)]l 

'  o  J 

(3.2 


where  we  have  used  the  fact  that  the  integrand  is  even  in  a  and  where 


X(n)  =1 


2  rZ 


Nc  -  P' 


1/2 


P2  -  f2 


(3.28a) 


<5"(n)> 


N 


<5‘(0)>  ^ 


(3.28b) 


For  completeness,  we  also  give  the  spectrum  in  (3,P)  space: 


AUTOSPECTRUM  ((CM)'/CPH) 


(c) 


FREQUENCY  (CPH) 


FIGURE  7 

TEMPERATURE  AND  UP  AUTOSPECTRA;  (a)  TEMPERATURE;  <b)  UP  (i.e., VERTICAL 
DISPLACEMENT  BASED  ON  THE  MEASURED  MEAN  VERTICAL  TEMPERATURE 
GRADIENT);  (c)  WKBJ  NORMALIZED  UP  AUTOSPECTRA. 

DIFFERENT  SYMBOLS  DENOTE  DIFFERENT  DEPTHS: 

SQUARE:  604m  CIRCLE:  731  m  TRIANGLE:  1023  m  CROSS.  2050  m 
(FIGURE  TAKEN  FROM  REF  34) 


ARL  UT 
AS-84-824 
CMD -GA 
9  ■  24  •  84 


39 


Chapter  4 


In  this  chapter  we  review  normal  mode  theory  for  a  horizon¬ 
tally  stratified,  deterministic  ocean.  In  particular,  we  consider  a 
single  water  layer  of  depth  H  overlying  a  fluid  half-space  of  constant 
compressional  velocity  Cp,  which  will  be  larger  than  the  sound  speed 
minimum  c  .  In  the  following  analysis  we  consider  a  time-harmonic 
source.  It  is  therefore  more  convenient  to  describe  sound  propagation 
by  means  of  a  velocity  potential  if (r » z )  such  that  the  particle 
velocity  u.  due  to  the  acoustic  wave  is 


u  =  V\|>(r,z)  e" 


(4.1) 


In  this  formalism,  the  acoustic  pressure  p  is  simply  related  to  by 


p  =  iwpil>  e 


(4.2) 


where  o  is  the  mass  density  of  the  medium.  Changes  in  density  with 
depth  are  closely  related  to  changes  in  the  sound  speed,  as  we  have 
seen.  Over  the  depth  of  the  ocean,  the  sound  speed  and  density  change 
very  little  compared  with  changes  in  >P .  Thus,  we  may  write  (from 
Eq.  (2.22))  a  wave  equation  for  . 


vS  +  ^  *  =  q(r,z) 
cc 


(4.3) 
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where  Q(r,z)^-^— represents  the  source  terms  which  are  the  right  side 
of  Eq.  (2.22).  The  simple  source  we  shall  consider  is  the  time- 
harmonic  point  source 


Q(r,z)  =  -4'rrQ05(^-x0) 


(4.4) 


where  Xq  is  the  position  of  the  source.  We  shall  choose  our  coordi¬ 
nate  system  such  that  xo=(0,zQ);  the  water-air  interface  is  at  z=0. 

In  the  deterministic  medium  chosen  here,  c  varies  only  with 


depth  so  that 


\p(r,z)  =  <J>(z)F(r) 


(4.5) 


That  is,  there  is  no  azimuthal  dependence  and  ^  separates  into  a 
product  of  functions  which  satisfy  ordinary  deferential  equations. 


In  particular 


(4.6) 


where  k  is  the  separation  constant. 

The  boundary  condition  that  the  normal  component  of  stress 
be  continuous  across  an  interface  implies  that  o<j>  is  continuous  across 
such  an  interface  and,  since  uz  must  also  be  continuous  across  an 
interface,  so  is  Since  the  mass  density  of  air  is  much  smaller 
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than  that  of  water,  we  may  set  <j>(0)=0.  Finally,  only  downgoing  waves 
are  permitted  as  z-*». 

If  k  is  greater  than  ■— ,  we  have  a  family  of  mode  functions 

which  exponentially  decay  as  z-+oo.  These  modes  correspond  to  waveguide 

phenomena  in  which  energy  is  trapped  by  total  reflection.^2  For  k>  — , 
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the  eigenvalue  spectrum  is  discretely  indexed  and,  if  there  are  M 
discrete  eigenvalues. 


—  >  kj  >  k2  >  k.  >  . 


k  >  - 

’  M  Cr, 


The  eigenfunction  <f>  . (z)  corresponding  to  the  ith  discrete  eigen¬ 
value  will  have  i  zeroes  in  the  depth  interval  [0,H). 

For  k<—  ,  the  spectrum  of  k  is  continuous.  The  modes 
CP 

corresponding  to  this  continuous  spectrum  are  generally  neglected 

because  they  are  attenuated  much  more  rapidly  than  the  discrete  modes, 

due  to  both  the  geometric  spreading  and  the  volume  absorption  which  is 

generally  smaller  in  the  water  layer  than  in  any  underlying  layers. 

It  is  a  simple  matter  to  show  that  the  discrete  eigenfunctions  ^ 
44 

are  orthogonal.  We  have  also  normalized  them  so  that 


dzptb  4  =  , 

n  m  nm 


(4.8a; 


dzp$(k)$(k' )  =  <5(k-k' ) 


:4.3b; 


I 


dzp(J>n4>(  k) 


(4.8c) 


Turning  to  the  range  functions,  we  find  upon  imposing  a 
Sommerfeld  condition  that 


F(r)  «  H^(kr) 


(4.9) 


where  H^(kr)  is  a  zeroth-order  Hankel  function  of  the  first  kind. 
Thus 


^(r,z)  =  iTrQop(zo) 


X^n{z)^n(zo)Ho1)(knr) 

n=l 


(- 


dk<j»(z,k)4>(20,k)Hj1^(kr) 


(4.10) 


(For  more  detail  see  Ref.  45).  From  now  on  we  follow  established 
practice  in  dropping  the  continuous  contribution. 

The  normalization  (Eq.  4.8)  has  a  very  specific  physical 
meaning.  Suppose  we  have  an  arrangement  of  sources  such  that  only  one 
of  the  discrete  normal  modes,  say  the  nth  mode,  and  none  of  the  others 
is  excited.  Then,  the  acoustic  Poynting  vector  is  given  by4^ 
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i"  J  (pu’+up*)  =  * 


q^o2(z0)p<z)^(z0)  (v[*n(Z)H<’>(knr)])’ 


’(>n(^)H‘T)(k„r))«n(z)H<Z)(k„r)]  . 


(4.11) 


The  integral  of  j[  over  depth  will  give  the  amount  of 
intensity  passing  through  a  line  at  horizontal  position  r. 


/  dzi  *  ?  F  «Jo 


p2(zo)4>n(zo^  I  dzp(z)^^z^ 


(4.12) 


For  this  cylindrically  symmetric  case,  the  intensity  passing 
through  a  cylinder  of  radius  r  is  thus 


Ir  =  2Ti2q2p2(zo)^(z0)u)  I  dzp(z)02(z) 


(4.13) 


If  we  define  range  functions 


Bn(0  =  i^0P(z0)^n(z0)Hj1)(knr) 


(4.14) 


then  the  velocity  potential  due  to  a  point  source  is 


=  £  Bn(r)*n(z) 


(4.15) 
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and  the  acoustic  Poynting  vector  is 


1  =  ?  7Qopo(zo)u)20n(zo)p(z)<J)n(z) 


n  ,m 


rfn(z>+m(z) 


Bn(D 


dB  (r) 
m 

dr 


Bjr) 


(r) 


dBm 

m 

TF 


2t„(z) 


d<>(z) 

vm 

~Sz 


(r)Bm(r)  - 


B*  ( r )  B. 


(r) 


(4.16) 


The  first  term  on  the  right  hand  side  of  Eq.  (4.16)  is  the 
sum  of  the  intensities  of  the  individual  modes  had  each  been  excited 
by  itself;  this  term  may  therefore  be  called  a  "self”  energy  term. 

The  restricted  sum  represents  the  intermodal  interference  and  is  truly 
an  interaction  term.  Note  that  the  self-energy  term  travels  in  a 
purely  radial  direction;  the  intermodal  interference  is  thus  solely 
and  explicitly  responsible  for  the  bending  of  ray  paths  in  the  ocean. 


Chapter  5 


Although  the  spectral  density  of  the  sound  speed  fluctu¬ 
ations  depends  only  on  the  magnitude  of  the  horizontal  wave  number, 
and  not  on  its  direction  (Eq.  3.25),  there  will  certainly  be  random 
azimuthal  fluctuations  in  the  sound  speed.  These  fluctuations  will  in 
turn  manifest  themselves  through  azimuthal  variations  in  the  acoustic 
pressure. 

fi-  ft 

Previous  studies  in  acoustic  normal  mode  theory  have 
neglected  this  azimuthal  variation  in  the  acoustic  field.  Kohler  and 

Q 

Papanicolaou  did  not  specifically  consider  internal  waves  as  the 
perturbation  mechanism;  they  simply  assumed  a  priori  that  the  sound 

C 

speed  fluctuations  did  not  depend  on  azimuth.  Dozier  and  Tappert, 
who  used  the  theoretically  derived  description  of  ^(x,i)(Eq.  3.17), 
integrated  the  spectral  density  of  the  sound  speed  fluctuations  along 
the  direction  perpendicular  to  the  line  connecting  the  acoustic  point 
source  and  the  receiver.  Therefore,  the  function  they  employed  as  the 
spectrum  is  actually  a  spectral  density  of  wave  number  components  in 
the  direction  of  propagation.  By  then  approximating  the  function 
Wj(a,z)  in  Eq.  (3.17)  by  Wj(o^,z),  where  is  the  component  of  a  in 
the  direction  of  propagation,  all  vector  dependence  by  5  on  horizontal 
position  is  replaced  by  a  scalar  dependence  on  range.  Let  be  the 
component  of  a  perpendicular  to  ct  •  The  justification  for  the 
replacement  of  Wj(a,z)  by  W^a^z)  is  that  "this  is  a  good 
approximation  for  small  a and  for  large  (and  hence  large  a)  the 


spectrum  will  be  low  anyway."  The  key  words  here  are  "small"  and 
"large"  since  they  immediately  raise  the  question  of  comparison. 
Clearly  the  comparison  is  between  and  this  approximation  is 

good  for  a ^  "small"  compared  with  ct^,  so  we  would  expect  the 
approximation  to  be  invalid  near  the  peak  of  the  spectrum  for  0.1  and 
both  small  and  about  the  same  size.  Thus  large  errors  in  the 
description  of  the  acoustic  field  may  or  may  not  result. 

The  most  disturbing  effect  of  this  assumption  is  to  preclude 
scattering  of  energy  among  the  azimuthal  acoustic  modes;  the  energy  is 
constrained  to  remain  in  the  azimuthal ly  symmetric  mode.  It  is  the 
purpose  of  this  chapter  to  discuss  this  case  but  without  making  any 
assumptions  about  the  relative  sizes  of  and  a?.  That  is,  we  allow 
azimuthal  fluctuations  in  the  sound  speed  but  assume  that  azimuthal 
coupling  in  the  acoustic  field  may  be  neglected.  We  will  thus  employ 
the  continuum  description  of  the  internal  wave  field  discussed  in 
Chapter  3.  In  addition,  we  do  not  make  the  "random  phase  approxima¬ 
tion";  different  depth  modes  may  indeed  be  correlated.  We  do  not 

expect  drastic  departures  from  the  results  obtained  by  Dozier  and 
6 

Tappert  for  this  case  where  we  have  artificially  prohibited  azimuthal 
coupling.  Still,  this  simpler  case  is  useful  in  examining  many  of  the 
physical  phenomena,  which  may  be  somewhat  obscured  by  the  added 
complexity  of  the  azimuthal  coupling,  and  in  introducing  the  method  of 
analysis  which  will  also  be  employed  later  when  we  do  not  ignore  this 


A.  The  Coupled  Mode  Equations 

Consider  a  flat  slab  of  seawater  of  thickness  H  overlying  a 
substrate.  The  velocity  potential  f(r,z,t)  due  to  a  time-harmonic 
point  source  of  angular  frequency  satisfies  the  scalar  wave  equation 
(cylindrical  coordinates) 


1  _3_  /  \  .  82V  1  92f'  . 

r  3r  \  dr  J  ^2  c2  gt2 


-4ttQ06()(_-x0)  e"iwt  _ 


(5.1) 


Because  the  sound  speed  fluctuations  vary  with  time  and  range, 
variables  do  not  exactly  separate.  However,  sound  can  travel  a 
distance  equal  to  the  correlation  range  of  the  fluctuations  (about 
10  km)  in  travel  times  on  the  order  of  5-10  sec,  much  less  than  the 


shortest  period  of  the  internal  waves  N 


-1 


20  minutes.  The 


internal  wave  field  is  thus  considered  "frozen"  with  respect  to  the 

1 6c 


acoustic  field.  We  still  have  the  range  variations  but  since 


«  1 


(see  Eq.  3.1)  we  can  expand  T  in  terms  of  the  normal  modes  4>n(z)  of 
the  deterministic  sound  speed  profile  c(z): 


:z) 


dz 


c2(z) 


-  V2> 


kn4n‘z' 


(5.2) 


with  the  normalization  (4.8a)  and  where  the  relative  densities  of  air 
and  water  allow  us  to  take 

;n(0)  -  0  .  (5.3) 
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If  the  substrate  is  rigid, 


<4n(H) 

cfz 


=  0 


(5.4) 


If,  on  the  other  hand,  the  substrate  is  an  isovelocity  fluid  with 
sound  speed  c^, 

*n(H)  =  Wn  F  exp 
w 

where  ff  is  determined  by  the  normalization,  P  is  the  water  density 
n  w 

at  z=H  and  P$  is  the  substrate  density  at  z=H.  We  do  not  consider  any 
other  boundary  conditions  in  this  study.  Thus 

M  _iujt 

¥(r,z,t)  =  ip(r,z)  e'Vjt  =  S  B  (r)4>  (z)  - -  ,  (5.6) 

n=l  Jr 

where  the  cylindrical  spreading  is  made  explicit. 

By  expanding  our  velocity  potential  in  terms  of  discrete 
normal  modes,  we  neglect  any  continuous  spectrum  of  modes.  In  the 
r'iid  bottom  and  pressure  release  bottom  cases,  there  are  an  infinite 
number  of  discrete  modes.  Where  real  bottom  cases  are  concerned,  M  is 
finite  and  there  also  exists  the  continuous  spectrum  of  modes  without 
which  the  normal  mode  set  is  not  complete.  The  traditional  assumption 
that  this  continuous  spectrum  does  not  contribute  greatly  to  the 


(5.5) 
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acoustic  field  is  made  here.  Therefore,  given  that  we  know  B  (r)  and 
4>  ( z ) ,  Eq.  (5.6)  may  not  be  strictly  true  but  it  is  certainly  a  good 
approximation.  A  more  important  error  incurred  by  this  approximation 
occurs  in  the  calculation  of  B  (r )  since  mode  coupling  into  the 
continuous  spectrum,  not  considered  here,  may  be  an  important  loss 

O 

mechanism.  Kohler  and  Papanicolaou0  have  advanced  a  formalism  to 
account  for  this  coupling  and  an  application  of  this  formalism  to 
sound  speed  fluctuations  caused  by  internal  waves  would  be  interesting 
future  work. 

All  of  the  randomness  in  the  acoustic  field  can  now  be  found 
in  the  range  functions.  Since  |6cj«c. 


w2  J  /,  26c\ 

q '  ?  v  ■  t) 


Employing  the  orthonormality  of  the  depth  functions,  we  find 


d  B  (r)  /  2  1  , 

^+(k"  +  i7)6"(r) 


if  dZo(z)V2>  ^377  »J2>6m(r> 

m  Jn  c  (z) 


Now  we  need  to  choose  initial  conditions  for  the  B  (r).  If  there  were 
no  perturbations  due  to  internal  waves,  the  functions  B  (r)  would  be 
the  expression  given  by  Eq.  (4.14)  multiplied  by  Jr  (see  Eq.  5.6). 


The  farfield  form  of  the  Hankel  function  is  accurate  within  4%  for 
knr>5  and  the  accuracy  increases  rapidly  with  argument.  The 
deviations  from  this  functional  form  due  to  the  internal  waves  is 
definitely  a  long  range  effect;  the  correlation  range  of  the  internal 
waves  is  about  10  km.^’B’^  It  is  therefore  convenient  to  separate 
the  ocean  into  two  regions:  a  "nearfield"  and  a  "farfield".  The 
nearfield  will  consist  of  the  region  interior  to  a  cylindrical 
boundary  centered  at  the  source  with  radius  rQ  =  (see  Eq.  4.7). 

”2  —  I 

Thus,  even  if  kM  ~  10  m  ,  an  order  of  magnitude  smaller  than  any  of 
the  eigenvalues  we  consider  in  this  study,  rQ  is  still  only  about  10% 
of  the  correlation  range  of  the  fluctuations.  Within  this  region,  we 
assume  that  Bn(r)  is  given  by  Eq.  (4.14)  multiplied  by  y?. 

In  the  farfield  (r>rQ),  we  neglect  the  l/4r  term  since  this 
is  less  than  1%  of  k^.  Bp(r)  will  satisfy  Eq.  (5.8)  with  initial 


condition 


/o-  i  (k  r  -tt/4) 

Bn(ro}  =  VlT  V(zo)<J,n(zo)  e 


In  terms  of  forward  and  backward  traveling  waves 


(5.9) 


i  /j. 
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(5.10) 
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The  quantities  A*  and  A”  are  complex  random  functions  and  we  are 
free  to  impose  one  relation  between  them.  We  choose 


iknr 


^n 

dr 


-’V 


dr 


=  0 


(5.11) 


so  that  Eq.  (5.8)  becomes 
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(5.12) 
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We  now  make  the  traditional  "forward  scattering 
approximation"  by  assuming  that  the  rapidly  oscillating  terms  on  the 
right  are  statistically  unimportant.  Also  called  the  "parabolic 
approximation"  because  of  the  form  of  Eq.  (5.12),  this  assumption 
is  equivalent  to  considering  the  small-angle  scattering  much  more 
likely  than  large  vertical  deflections.  The  evidence  for  this  is 
mainly  exper imental  ^  any  ra^e#  the  s-jze  0f  1S 

generally  several  orders  of  magnitude  smaller  than  (kn+km).  It  will 
be  shown,  consistent  with  the  results  of  Dozier  and  Tappert^  and  of 

o 

Kohler  and  Papanicolaou,  that  the  effect  of  the  argument  in  the 


exponential  is  to  pick  out  from  the  internal  wave  field  scatterers 


whose  horizontal  vector  wave  numbers  have  a  component  the  size  of  the 
arguments  of  the  exponentials  in  Eq.  (5.12).  Since  the  internal  wave 
spectrum  is  heavily  weighted  toward  small  horizontal  wave  numbers, 
it  is  much  easier  to  find  scatterers  with  horizontal  wave  numbers  on 
the  order  of  (kn-k  )  than  those  with  horizontal  wave  number  (kn+km). 
We  therefore  ignore  the  backscattered  wave  and  (dropping  the  plus 
sign,  consistent  with  A"  being  initially  zero), 


~  =  icY  V  A 
r  / .  nm  m 


(5.13) 


where 


7^=r/  dz‘>(z)<'n(2>  »m<z>  e 

eVknkmJo  c 


i  (k  -k  )r 

m  n  .  (5.14) 


e 

The  quantity  An  is  the  envelope  of  the  rapidly  oscillating  —  term 
and  would  be  constant  in  a  deterministic  medium.  Since  we  estimate 
that  changes  of  order  unity  will  occur  in  Ap  over  a  range  of  about  a 
correlation  range  of  the  fluctuations  (C^  ~  10  km)  we  have  introduced 
the  smallness  parameter  £  in  order  to  make  Vnm  of  the  same  order  as 
1/CR.  Our  task  now  lies  in  the  estimation  of  e. 

From  Eq.  (5.14),  if  Vnm~l/CR,  k-k^u./^,  and  c-cQ,  then 


-af 

EcoJ„ 


dzp(z)Cm(z)  —  4n(z) 


(5.15) 


,  and  also 


2  1/2 

If  we  now  replace  by  the  depth  average  of  <(-=■)  > 
assume  that  the  mode  functions  contribute  very  little  for  z?H, 


-?(*)<  (»t, 


1/2 


(5.16) 
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ft  //J  p  \  *“  *  I  *- 

Dozier  and  Tappert  estimate  <(^£-)  >  ~  5.6x10"  ,  so 

'c  7z=o 


c  ~  (4.04  xio-4) 


w 


(5.17) 


The  pertubation  theory  is  thus  seen  to  be  a  low  frequency  theory 

and  should  break  down  for  greater  than  a  couple  of  hundred  Hertz. 

The  expression  (5.16)  for  e  reveals  that  our  smallness 

parameter  is  really  a  product  of  two  competing  considerations,  the 

size  of  the  sound  speed  fluctuations  and  also  their  coherence.  The 
/  B  \  /<Sc\2  1/2 

quantity  j  is  an  estimate  of  how  big  the  fluctuations 

are  compared  with  the  deterministic  sound  speed.  Clearly,  if  a 

sizable  fraction  of  the  sound  speed  is  due  to  the  fluctuations,  then 

U)C 

Eq.  (5.7)  will  not  be  valid.  On  the  other  hand,  the  quantity  — - 

co 

compares  the  correlation  range  of  the  sound  speed  fluctuations  with 
the  acoustic  wavelength.  If,  in  the  horizontal  plane,  the 
characteristic  length  of  the  sound  speed  fluctuations  is  only  a 
couple  of  acoustic  wavelengths  long,  the  sound  speed  fluctuations  are 
not  strongly  correlated,  and  the  acoustic  wave  will  not  travel  very 
far  without  being  effectively  sped  up  and  slowed  down  several  times  by 
the  random  effects  on  the  phase.  As  the  phases  get  more  and  more 


mixed  up,  the  acoustic  wave  will,  in  a  sense,  forget  where  it  has 


been.  On  the  other  hand,  a  large  CR  means  that  internal  wave  effects 
on  acoustic  propagation  will  be  more  systematic.  These  effects  will 
then  tend  to  accumulate  and  knowledge  of  past  history  will  be  more 
important  than  in  a  system  where  the  sound  speed  fluctuations  are  the 
same  size  but  less  correlated.  Since  the  correlation  range  is  many 
times  the  acoustic  wavelength,  this  pertubation  theory  is  valid  only 
because  the  sound  speed  fluctuations  are  so  small. 

B.  The  Statistical  Behavior  of  the  Range  Functions 
1.  Interesting  Moments 

The  randomness  in  the  sound  speed  translates  into 
randomness  in  the  velocity  potential  and  any  function  f(if(r,z))  of  the 
velocity  potential.  We  are  interested  in  moments  of  these  random 
functions. 

The  quantity  <  (r,z) >  is  definitely  not  the  value  of 

4>  in  the  absence  of  any  perturbation.  If  Nature  had  simply  added  a 
random  term  to  the  wave  equation  rather  than  multiplied  it  by  the 
state  of  the  system,  then  first  moments  of  quantities  like  the 
acoustic  pressure  would  indeed  be  their  deterministic  values,  although 
higher  moments  would  be  somewhat  more  complicated.  However,  the 
nat  're  of  our  problem  is  explicitly  multiplicative  and  it  is  thus  of 
some  interest  to  calculate  the  first  moment  of  ^(r,z). 


i 


J 


K 


In  the  perturbative  method  used  here,  we  will  calculate 
moments  up  to  lowest  non-vanishing  order  in  e.  That  is,  any  average 
quantity  can  be  ordered  with  e : 


<f 


(Mr,z))>  =  fQ  +  t  <f]  >  +  c2  <f2>  +  e3<f3>  +  ...  ,  (5.18) 


where  f  is  the  initial  value  of  f.  If,  for  example,  c  <f^> 
vanishes,  we  wi  1 1  calculate  express  ions  for  <f>  accurate  to  O(e^). 

Three  quantities  in  which  we  are  particularly  interested  are 
the  first  moment  of  the  pressure 


<D>  .  l,,0£  <An(r)>  »n(z> 
n=l 


P>  =  iiopj - - e  "  ,  (5.19) 


A7 


the  second  moment  of  the  pressure 


.  .,,2^  _  2  2V  V  »n<zH„<z>  J<Vkm)r 

<|p(r,z,t)  |  >  =  tu  p  y  y  - - - - e 

n=1  r  AS 


and  the  first  moment  of  the  acoustic  Poynting  vector 


(5.20) 
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(5.21) 
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Equations  (5.19-5.21)  are  not  exact  expressions  for 
these  quantities  since  we  have  already  dropped  some  higher-order  terms 
in  using  Eq.  (5.7).  At  any  rate,  we  shall  content  ourselves  with  the 
lowest-order  expressions  for  <p>,  < [ p [  >,  and  <_j>.  Thus,  we 
drop  the  term  proportional  to  e  in  Eq.  (5.21).  In  addition,  since  we 
are  considering  farfield  propagation,  we  also  drop  the  l/2r  term. 
Hence, 
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(5.22) 


Note  that  again  <j>  is  divided  into  self  and 
interaction  energy  terms.  The  autocorrelation  of  a  mode  with  itself 
represents  the  self-energy  while  interference  effects  are  represented 
by  the  energy  contained  in  correlations  between  the  modes. 

2.  Calculation  of  the  Moments 

The  scaling  method  used  here  to  estimate  the  moments  of 
A,  where  A  ,  was  introduced  by  Stratonovich.  ^ 


The 


literature  extending  this  method  to  increasingly  larger  classes  of 
problems  is  vast  (see  Ref.  50  and  references  within);  the  version  of 
interest  here  was  advanced  by  Papanicolaou  and  Kohler.^  Although  the 
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work  by  Kohler  and  Papanicolaou  and  Kohler  seem  to  be  the  only 
studies  of  the  effects  of  internal  waves  on  acoustic  normal  modes  to 
employ  rigorous  scaling  methods,  these  methods  have  been  used  in  com¬ 
munications  research  for  many  years  to  describe  wave  propagation 
through  waveguides  with  random  inhomogeneities  (e.g.,  Refs.  52  and 
53). 

The  method  entails  scaling  the  range  coordinate  by  the  size 
of  the  fluctuations.  In  this  scaled  coordinate  ?*£  r,  the  effects  of 
the  random  fluctuations  are  squeezed  together  enough  that  the  acoustic 
field  can  be  approximated  by  a  Markov  diffusion  process  (Appendix  C) 
with  an  error  of  0(e'  in  the  scaled  coordinates,  or  O(s^)  in  the 
original  coordinate  system.  There  are  two  great  advantages  to  this. 
First,  operational  methods  of  treating  Markov  diffusion  processes  are 
available.  Second,  we  have  an  estimate  for  the  error  in  terms  of  its 
order  in  c. 

This  approximation  will  hold  at  long  ranges  for  small  values 
of  - .  Thus  from  Eq.  (5.13) 

dA 

d ,  (5.23) 

where  _V=  •  As  e~0,  r-^  witn  :=-cr  fixed,  moments  of  A  converge  to 


We  assume 


correspond ing  moments  of  a  Markov  diffusion  process  A0, 
that  the  sound  speed  fluctuations  at  two  different  ranges  in  the 
unsealed  coordinate  system  become  increasingly  independent  in  a 
sufficiently  strong  sense  (for  a  definition  of  "sufficiently  strong" 
see  Papanicolaou  and  Kohler^)  as  the  range  difference  increases.  We 
also  demand  that  ^  ^  finite  but  arbitrarily  large.  The 

limitation  of  a  finite  has  been  1  i f ted  for  stationary  systems."'' 
However,  our  system  is  not  stationary  in  range,  but  in  horizontal 
position.  Thus,  we  are  confined  to  the  weaker  results  of  Papanicolaou 
and  Kohler.^  This  does  not  present  a  problem  since  the  ocean  isn't 
really  infinite  anyway. 

The  system  (5.23)  is  completely  analogous  with  the 

10  ? 

example  in  Remark  2  of  Papanicolaou  and  Kohler;  if  c0=c  r0> 
v(A°,  A°  )  is  a  function  of  the  limiting  Markov  process  A0, A0*,  and 


$4  ;~o!  v(«°,A°  ) 


(5.24) 


where  p(A°,A°  ,A^;cq)  is  the  transition  probability  density  of 

finding  the  system  in  state  A°A°  at  scaled  range  "  given  initial 
conditions  A  ,A°  (aQ)=A°,A^*  then  T  (A°,Aq*)  obeys  a  backward 
eguation, 


where 


t  =  Tin 
T** 


nip  d,j  dt  <I(1(,)S#'^)*(,I 
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(5.261 


In  Eq.  (5.26),  ( — ^ \  =  — rr  and  the  inner  product  (A,B)  is  defined  as 

\ml  ^Ano 


(A ,3)  =  tr(AB  )  . 


(5.27) 


Since  limits  have  been  placed  on  the  amplitude  of  the 
acoustic  signal  (or  else  we  can't  start  out  with  the  linearized  wave 
equation),  we  assume  that  '(A0, A0*)  vanishes  for  A°  ,A°*  larger  than 
some  bounds.  If  we  also  assume  that  first  and  second  derivatives  with 
respect  to  modal  amplitude  exist,  then  the  transition  probability 

(0  0*  i  0  0*  \ 

A  ,A  ;c|A 0,A0;aQj  obeys  a  forward,  or  Fokker-Planck , 

equation 


dp(A°,A°VAX*;, 


=  t  p(A0,A°*;a!Aj,Aj*;c0  )  ,  (5.23) 


where  if,  is  the  formal  adjoint  of  the  operator  defined  in  Eq.  (5.26) 
and  p(A°,A°*;r0'Aj,Aj*  =  c^A°-^j :  ^A°*-A°*  j.  One  additional 
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caveat:  in  order  for  the  scaling  approximation  to  hold,^  must  exist 
independently  of  a  .  Papanicolaou  and  Kohler^  assumed  this  in  order 
to  relax  the  condition  that  the  fluctuating  medium  be  stationary  in 
range.  There  is  no  real  way  to  know  whether  is  really  independent 
of  aQ  for  a  nonstationary  process  without  actually  calculating  it. 
This  is  not  a  problem  when  we  neglect  azimuthal  coupling  in  the 
acoustic  field,  but  we  will  need  to  remember  this  later  when 
azimuthal  coupling  is  taken  into  account. 

Equation  (5.28)  is  very  useful  in  calculating  propagation 
equations  for  the  moments  we  need.  Since  we  are  approxi¬ 
mating  our  real  process  by  the  limiting  process,  we  will  henceforth 

drop  the  superscript  on  A0.  By  multiplying  Eq.  (5.28)  by  A^  and 

★ 

integrating  over  A  and  A  ,  we  get  a  propagation  equation  for  <  Ak  > 
(Appendix  0) : 
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Hence,  since  e  r  =  o,  we  find 


,  -Z  (amk-ib  ,)(r-r  ) 

<Ak(r)>  =  e  m  0  <Ak(ro)>  .  (5.32) 


The  details  of  this  calculation  are  given  in  Appendi/.  D. 

In  a  similar  manner  we  derive  propagation  equations  for 
the  second  moments: 
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so  that 
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.fca.i.'-.l.  1  .r 
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Even  with  the  oscillatory  terms,  a  .  >0.  Thus  <Ak(r)>  and 
<AkAk,(r)>  oscillate  with  an  exponentially  decreasing 

2 

envelope.  The  quantities  <|Akr  >  apDroach  an  equilibrium  value 


<iAki2>  *fil  <ifl™<ro>i2> 


(5.36) 


Hence,  knowledge  of  initial  conditions  is  eventually  wiped  out  as  the 
energy  is  equipartitioned  among  the  modes.  From  Eq.  (5.9) 
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so  that  the  equilibrium  values  for  <p>,  <|p(r,z,t)|  >,  and  <j> 
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Chapter  6 


In  Chapter  5,  we  studied  sound  propagation  where  each  depth 
mode  was  considered  to  have  an  azimuthal ly  symmetric  amplitude  for 
each  realization.  This  directly  contradicts  the  assumption  of  azi¬ 
muthal  variations  in  the  sound  speed  profile.  In  this  chapter,  we 
present  a  formalism  for  dealing  with  the  azimuthal  variations  in  the 
sound  field. 

A.  The  Coupled  Mode  Equations 

We  consider  the  same  system  as  in  the  previous  chapter,  a 
flat  slab  of  seawater  of  thickness  H  overlying  a  substrate.  The  veloc¬ 
ity  potential  T(r,ip,z,t)  satisfies  the  scalar  wave  equation 
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The  sound  speed  profile  is  again  considered  "frozen"  and  we  again 
expand  7  in  terms  of  the  deterministic  depth  modes  <j>n(z): 
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We  assume  that  internal  wave  effects  have  negligibly 
disturbed  the  acoustic  field  interior  to  a  cylindrical  boundary  at 
ro=10/kyj.  Thus,  Eq.  (6-3)  will  describe  the  exterior  problem  with  Bn(rQ) 
obtained  from  Eq.  (4.14)  multiplied  by  ,/r^.  Accordingly,  we  drop 
the  l/4r^  and  make  the  forward  scattering  approximation 

Bn(r,<p)  =  -4=-  B^(r,<p)  e  n  ,  (6.4) 
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Now  we  expand  the  forward  scattered  wave  in  terms  of  azimuthal  modes 
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B.  Calculation  of  Moments 

As  in  the  last  chapter  we  will  be  interested  in  expressions 
for  the  average  pressure. 
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the  second  moment  of  the  pressure. 
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and  the  average  energy  flux, 


We  will  find  that,  as  expected,  a  cylindrically  symmetric  source  will 
require  q=0  in  Eq.  (6.14)  and  though  q  need  not  be  zero  in  Eqs.  (6.15) 
and  (6.16)  it  must  be  true  that  q=p. 

Equation  (6.9)  for  the  matrix  process  A=(AnC|)  is  completely 
analogous  to  Eq.  (5.23)  for  the  vector  process  considered  in 
Chapter  5.  We  shall  place  a  restriction  that  q  be  less  than  some 
limit  L  where  L  is  arbitrarily  large.  That  is,  we  assume  that  we 
can  always  approximate  a  delta  function  by  a  function  that  is  highly 
peaked  at  <P=0,  the  width  of  which  is  small  compared  with  measured 
azimuthal  variations  in  the  velocity  potential,  and  that  we  obtain 
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this  function  simply  by  adding  together  an  adequate,  but  finite, 
number  of  azimuthal  modes. 


Given  this  restriction,  we  proceed  in  a  manner  analogous  to 
the  situation  in  Chapter  5.  Rewriting  Eq.  (6.9), 

dA 

H?  =  M  .  (6-17! 

where  ^=(Anq)  and  the  sca^n9  a=£^>  the  processes 

A(~y)can  be  approximated  by  the  Markov  diffusion  processes,  A °(o): 


dp(4°,A°ta|^0,A^*a0) 


da 


■=£p(A°,A°*;a 


(6.18) 


whereat'  is  the  formal  adjoint  of 


The  form  of  V  (Eq.  6.10)  demands  that  our  ao  be  nonzero  (which  it  is) 

-L  A. 

in  order  to  calculated  but  we  will  finddTto  be  independent  of  aQ 
when  we're  finished.  In  terms  of  the  components  of  V  and  A  (dropping 
the  superscripts), 


,  ,  ra.+T  rS 

£  =  2  lim  , 
fppi?q ,  ,  Tf0°  '  JoQ 
ni  p  n  q  u 


}JVT«f 


dt 
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6.6  , 


9  9  *  * 

t - V-  a__a. 


mpnq  m  p  n  q  \  n'mpq'  -  ^  -  "nq"n'q' 

m  p  m  p  mD 


*  ^V(t)  V*(s)  -  3  3  a*  a 

mpnq  m' p ' n 1  q ‘  >  —  An'qlAnq 

m'p'  dmp 


+  ,  ,>  — jr- 

mpnq  m'p'n  q  9ft 


m'p'  9A* 

K  '  mp 


(6.20) 


Each  of  the  four  terms  will  contain  an  expression  of  the  form 

fee  f  271  f2*  rrr/2  . 

J  dot  J  d<?  J  d(p' J  deP(a,n; n)  cos(a^x)  elAq<;p  elAq,<p' 


^e^[ocos(«p-0)-tc°s(<P'-9)]  +  eia[ocos(<A+6)-tcos(tp'+e)]J 


(6.21) 


By  transforming  <P-*-cP+0  ,  c^uco'  +  Oin  the  first  integral  and  <p-»(p-0  , 
<p' -KP'-0in  the  second  integral,  the  expression  (6.21)  becomes 

f"  r  n/2  r2TT  f2n 

J ^  daj  <16  J  d vj  dVP(a,r.;n)  cosh?X)  eUo,:°s,-\-'iatcosv 


-  A 


'  '  1 
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^  giAq(<P+e)  eiAq'(<p'+e)  +  eiiq(«P-e)  ei:'.q' 

?  sin["^(Aq+Aq')]  /*“  r2- r  /*2tt 

- (7q+Aq 7 ) - dot  J  d(?J  d^’ P( j ,7;r>)  cos(^X) 


<  eiaocosy»  e-iatcos^*  gi  (Aq<p+Aq  V 


(6.22) 


If  Aq  and  Aq 1  are  both  even  or  both  odd  the  integral  is  zero.  Even 
more  is  required:  because  P(a,f.;n)  is  even  in  a  and  because  the  angle 
integrals  go  from  0  to  2tt,  the  expression  (6.22)  is  null  unless 
Aq=-Aq'  (Appendix  E).  Hence, 


-Cn+T  f  S 


lim  j  I  ds  I  dt 
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■*Vd'  ‘  m '  p '  '  mp 


+  <M^  V*(S1  ,  ,>  .  9  9  fl*  A 

mpnq  m  p  n  q  np ,q+p •  _q •  7^—  An'q'Anq 
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+  <V*(t)V(^t  >  5  -  - _ iL_  a  A* 

mpnq  m  p  n  q  ^p.q+p'-q'  3A  ,  ,  ;A*  Hn'q’  nq 

Mmp 


.  (6.23) 


Another  property  common  to  the  four  terms  in  Eq.  (6.18)  is  an 
expression  of  the  form 


rdaf°0+Tds  f  SdtP(o,:!in)  cos(a^X)  eMscow>-tcos«>' ) 

*'-00  Jrt  _ 


'°0  a0 


iAk.s  iAkpt  (i/2s )&Q^  (i/2t)AQ? 

x  e  e  e  e  ^ 


(6.24) 


where,  for  example  in  the  first  term,  Ak^  =  ( k  -k  ) ,  Ak^= ( km • -k  . ) 


/  2  2\  9  /  ,2  i 2 \ 

=(F  -  and  aq2  =  (E—  -  %—)  •  s 

\km  n  f  2  \V  V  / 


ince  -j— cos(a^  X)  <l/a*, 
oi  +a* 


AQ^  and  AQ^  are  constrained  to  be  finite  although  arbitrarily  large,  oQ 
is  finite  (not  zero)  and  less  than  T  but  otherwise  arbitrary;  and  since 
no  function  of  SI  multiplies  T  in  the  argument  of  an  exDonential,  the 
expression  in  (6.24)  can  be  written  (Appendix  F) 


i  H  da  f^ds  fS  dtP(a,Q;n)  COs(a^X)  eMscosp-tcos<P'  )  e1Akls  J' 


m 


Since  the  terms  o(lnT/T)  and  0 ( 1/T )  will  vanish  in  the  limit  T-*»,  they 
are  henceforth  omitted. 

It  is  now  possible  to  calculate  some  useful  moments.  In  the 
same  manner  as  in  the  last  chapter,  a  propagation  equation  for  <  j > 
is  obtained  by  multiplying  Eq.  (6.18)  by  and  integrating  by  parts 
once.  Thus 


d<\j> 
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(6.26) 
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Now 


e1p(«’-'1’,)  *  2M to.-®') 

P=-L 


(6.27) 


where  this  is  approximate  because  of  the  truncated  sum.  Further 
manipulation  yields 


d  <Ak.i > 

do 
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Jo  Jq  Jo 


dz'^(z)^m(zHm(z'Hk(z') 


P(Z)0(Z') 


/•  M  ( n )  po 
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3~ 


da_*  isl£i 

(a  +a$)  sr 


x  f  °  ds  [S  dt  el(km'kn)s  el(kk‘km)t  eiacosa,(t-s)  <Ap.>  . 

(6.28) 

Notice  that  the  operator  on  <Anj>  does  not  depend  on  the  indices  of 
any  of  the  azimuthal  modes.  In  fact,  this  operator  is  precisely  the 
one  evaluated  in  Appendix  D  and 


d  <Ai.i > 


<Vi> 

E  m 


(6.29) 
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with  a^  and  t>mk  given  by  Eqs.  (5.30)  and  (5.31).  Thus 


xL  (amk'ibmk) (r-ro^ 


16.30) 


The  derivation  of  Eq.  (6.29)  does  not  involve  the  initial 
conditions;  the  operator^  is  purely  a  property  of  the  medium  and  the 
acoustic  frequency.  Hence,  Eq.  (6.30)  holds  for  directional  sources 
as  well  as  azimuthally  symmetric  sources;  the  first  moment  preserves 
the  original  directionality  of  the  source.  This  is  expected  since  the 
first  moment  does  not  give  any  information  about  energy  transfer  but 
merely  serves  as  an  indication  of  how  quickly  knowledge  of  the  initial 
conditions  is  wiped  out. 

Information  about  the  spatial  distribution  of  energy  is 

★ 

obtained  from  the  second  moments  <  AnqAn  <  q ■ >  (Appendix  G).  Here,  it 
is  necessary  to  return  to  the  azimuthally  symmetric  source  in  order  to 
make  progress.  We  find  first  of  all  that  the  different  azimuthal 
modes  are  uncorrelated  to  lowest  order  in  £ ;  that  is,  we  must  have 
q=q'.  The  random  pressure  p(r,<P,z,t)  is  thus  stationary  in  the  hori¬ 
zontal  plane. 

The  interesting  quantities  involving  second  moments  are 


<l  ,2,  V  C(z)C’1;)  ei(kn'kn',r 

r  AaT 


<  A  A  ,  >  6X  D 
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,  (6.31) 
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<  j  > 
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,,  2  2 
2kn,r 


+  3.  -Jl——--  Re  <A  A  >  exp  ji  (kn-kn,)r  +  ^  ^ 

Aikn'  L  LL  n  n 


</Vn'q>  exp  VL(kn-kn,)r  +  §F  (k^ 

'v'  nKn'  L 


(6.32) 


It  is  necessary  here  to  make  some  estimate  of  our  limit  on 
q.  The  sound  speed  fluctuations  have  a  correlation  distance  of  about 
10  km.  Therefore,  traversing  an  arc  centered  at  the  source,  we  expect 
changes  of  order  one  in  the  acoustic  field  over  an  arclength  of  about 
a  correlation  distance.  We  thus  take  L  =  ( 2~r r ) / ;  our  limit 
increases  with  range.  Considering  the  radial  term  in  <j>. 


J.  r 
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<  2*10" 


(6.33) 


Simi 1 ar  ly 


L2  /  1 


2r  k 


-  rM;  M0-3(kn,-k„), 


n  n 


(6.34) 


Hence,  we  'nay  drop  the  V~ T  term  in  <  j  >  and  can  probably  also  drop 

2k  r  ~ 

n  2 

the  extra  oscillatory  terms  (6.34)  in  <|p|  >  and  <j_>,  depending  on 
the  size  of  (kni-kn)r.  In  the  estimates  given  above,  k  ~10  2  rr  ^ ,  about 
an  order  of  magnitude  smaller  than  what  k  would  be  at  an  acoustic 
frequency  of  100  Hz. 

★ 

The  quantities  ^  A  An  ,  ^  >  satisfy  the  following  equation: 
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(6.35) 


where  amn  and  bmn  are  given  by  Eqs.  (5.30)  and  (5.31).  The  sums  over 


m  and  m'  are  restricted 


such  that  !kn-km).<)knI-km( )  and  that  (k^kj 


have  the  same  sign  as  (k  The  quantities  dpqm,n,  are  given  by 
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(6.36) 
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Note  the  similarity  between  d , n ,  and  amn  (Eq.  5.30).  The 

excitation  of  the  azimuthal  modes  by  the  internal  wave  field,  bleeding 

energy  out  of  the  azimuthally  symmetric  mode,  is  explicit  in 

Eq.  (6.35).  Since  all  the  energy  is  initially  in  the  q=0  mode,  we  see 

from  the  form  of  d^,  that  <AnqAn,q>  =  <  An>_qA* ,  ^  >  at  all 

ranges.  Therefore,  <  j  >  does  not  have  a  component  in  the  Ip  direction. 

The  moments  <AnqAn,q>  are  too  many  and  too  complicated  to 

evaluate  numerically.  However,  we  obtain  much  information  about 

2 

energy  flow  without  calculating  each  <Ap  An,  >.  Both  <  |  pj  >  and 


<j>  contain  terms 


(6.37) 


which  obey  a  propagation  equation, 


=  2 }  a  fw  -w  1 
'  -  rrnv  m  n ' 


(6.38) 


2 

This  is  precisely  the  same  equation  as  that  satisfied  by  <  | A  |  > 

(Eq.  5.35)  for  the  case  when  azimuthal  coupling  was  neglected.  The 
self-energy  is  not  affected  by  modal  interference  and  so  is  simply  the 
sum  of  the  self-energies  contained  in  all  the  azimuthal  modes. 

The  interference  terms  are  another  story.  For  n^n’, 
summation  over  q  in  Eq.  (6.35)  yields 
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>]I 


<  AnqAn'q> 
(6.39} 


which  is  analogous  to  Eq.  (5.33)  for  the  azimuthal ly  symmetric  case. 

2 

However,  the  expressions  for  <|p|  >  and  <j>  contain  expressions 

which  oscillate  more  quickly  than  £  <-AnaAn'a>: 

q  q  q 
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<flnqfln'q>  exp 


la! 

2  r 


/ 1  1  \ 

< 

(  k  k  ,  1 

\  n  n  / 
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<^nq^n  'q> 


(6.40) 


The  effect  of  azimuthal  coupling,  then,  enhances  the  decorrelation  of 
the  depth  modes  so  that  the  interaction  terms  go  more  quickly  to  zero. 

This  enhancement  is  probably  not  very  big.  Changes  in  phase 
are  much  more  rapid  in  range  than  in  azimuth  and  a  slight  advancement 
or  retardation  of  a  depth  mode  as  it  travels  in  range  will  play  much 
more  havoc  with  acoustic  interference  than  lateral  deflections  of 
energy.  Still,  the  consideration  of  the  extra  dimension  can  only 
serve  to  increase  randomness  in  the  phases;  this  increase  is  repre¬ 
sented  by  the  additional  oscillations  in  the  correlations  between 
different  depth  modes. 

In  view  of  the  above  considerations,  the  average  energy  flux 
is  taken  to  be 
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Chapter  7 


In  order  to  determine  the  rate  of  acoustic  energy 

redistribution  by  the  internal  waves  and  the  decay  rate  of  the 

correlations  between  the  different  mode  functions  it  is  desirable  to 

numerically  calculate  the  scattering  coefficients  ann>-  Because  of 

the  inequality  (6.40),  it  is  certainly  necessary  to  know  the  decay 

rate  of  2_,  <■  A  A  ,  >  since»  if  this  quantity  be  negligible,  the 
9  2 

interference  terms  in  <|p|  >  and  <j>  may  be  neglected. 

The  scattering  coefficients  a  .  are  given  by 
a  nn 13  J 
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f2)V2 


and 


(7.5) 


(7.6) 


The  -  integration  can  be  transformed  by  a  change  of 


variable. 
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where 


Ak> 


(7.8) 


coscp 


K  =  N(n)/f 


(7.9) 


,2  _  t2(N2-f2)  cos2<P 
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Ak 


For  Ak>ft,  we  can  separate  the  denominator  of  the  integrand 
in  (7.7)  into  partial  fractions  and  evaluate  the  X  integral  exactly. 
For  Ak<ft,  the  partia1  fraction  expansion  will  introduce  second-order 


poles  at  ft  cos  </>=  Ak  and  then  the  principal  value  of  the  cp  integral 
will  not  exist.  It  is  certainly  true  that  the  expression  (7.1)  for 
ann,  is  valid  uniformly  in  latitude  and  bandwidth  parameter  t,  but  one 
must  be  careful  that  methods  chosen  to  evaluate  the  integrals  be 


valid. 


Evaluation  of  ann,  is  greatly  simplified  by  assuming  Ak>ft. 


This  limits  the  latitudes  at  which  the  following  expressions  apply, 

but  this  does  not  prevent  a  useful  computation.  For  t  =  1.4  sec/m, 
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consistent  with  j*=3  in  Desaubies'  0  formalism,  and  a  latitude  of  30°, 
-4 

ft  =  1.02x10  rad/m.  At  any  rate,  the  expression  (7.1)  is  as  refined 
as  it  can  get  without  introducing  some  approximations.  Numerical 
evaluation  of  the  quadruple  integral  for  each  pair  of  depth  modes 
would  take  a  prohibitive  amount  of  computer  time,  even  with  Monte- 
Carlo  techniques  of  integration.  Since  we  do  not  expect  the 
qualitative  behavior  of  the  acoustic  field  to  vary  drastically  with 
latitude,  we  solve  the  problem  for  Ak>ft. 

Performing  the  internal  wave  frequency  integral,  we  find 
(with  N»f;  Appendix  H) 
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(7.11 


where  Ei(x)  is  the  exponential  integral. 

Because  the  exponential  integral  diverges  near  ^  =0,  it  is 
interesting  to  note  that  as  Gr*-Q  the  terms  in  curly  brackets  become 


(7.12 


a  result  which  can  also  be  obtained  by  integrating  the  integral  (7.7) 
with  G=0  (Appendix  H). 

The  expression  (7.11.)  for  ann,  looks  much  more  daunting  than 
it  actually  is.  The  integrand  is  actually  very  smooth  and  well- 
behaved  except  near  ^  =0.  We  can  eliminate  much  repetition  in  the 
numerical  evaluation  of  a  ,  by  noting  that  the  expression  (7.11)  is 
of  the  form 

fH  fH 

ann,  =1  dz  dz'h(z)h(z' )g( f  ,n)  ,  (7.1 

*0  () 
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which  is  equivalent  to 


a 
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After  one  more  transformation 
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(7.15) 


The  scattering  coefficients  ann,  were  calculated  for  n=l, 
n ' =2-6  at  acoustic  frequencies  of  100  Hz  and  50  Hz.  We  assume  an 
exponentially  decreasing  Brunt-Vai sala  frequency 


N(z)  =  Nq  e~z/B 


(7.16) 


where  Nq  =  0.00524  rad/sec  and  B  =  1300  m.  Correspondingly,  a  Munk 
profile  was  chosen  for  the  deterministic  sound  speed 

c(z)  =  cQ[l  +  e(A+e'A-l)]  ,  (7-17) 

where  cQ  =  1500  m/sec,  e=0. 737x10’^, 

A=|(z'zo)  5  (7.18) 

-4 

and  zQ  =  1300  m.  The  inertial  frequency  was  taken  to  be  0.727x10 
rad/sec  (30°  latitude).  The  bandwidth  parameter  t  was  chosen  to  be 


■DC 

1.4  sec/m,  consistent  with  j*=3  in  Desaubies1  notation.  Finally, 
in  order  to  compare  our  numerical  results  with  those  of  Dozier  and 
Tappert,7  we  choose  the  strength  of  the  fluctuations  such  that 

(|)2  <r2{ o)>  :i0  -  2.12  s3  .  (7.19) 

The  acoustic  wave  numbers  kn  and  the  corresponding  depth 

modes  <t>  were  calculated  for  a  rigid  bottom  using  the  normal  mode 
54 

model  NEMESIS  on  a  CYBER  171  computer.  Briefly,  boundary  conditions 

were  set  at  z=0  and  z  =  H  =  4000  m.  Given  an  eigenvalue  estimate, 

Numerov's  method  is  used  to  integrate  the  differential  equation  (5.2) 

up  from  z=H  to  a  match  point  and  down  to  that  match  point  from  the 

surface.  The  difference  in  these  two  solutions  at  the  match  point,  if 

not  within  the  required  tolerance,  is  then  used  to  determine  the  next 

eigenvalue  estimate.  Polynomial  expressions  involving  previously 

calculated  eigenvalues  are  used  as  estimates  for  the  next  eigenvalue 

so  that  convergence  is  usually  achieved  in  five  or  six  iterations. 

Each  eigenfunction  <p  was  evaluated  every  5  m  so  that  there  were  801 

mesh  points  in  the  4000  m  of  water. 

Because  the  <P  integrand  is  so  smooth  and  since  the 

<p  integral  must  be  evaluated  for  each  (n  ,  O  )  pair,  the  integral  over 

azimuth  was  found  to  be  evaluated  most  efficiently  using  a  2*  point 

55 

Gaussian  tensor  product  formula,  where  I  varies  from  1  to  8  until 


55 

the  •-■eti'j  i t  ed  accuracy  is  achieved.  The  IMSL  routine  DMLIN  was  i  sed  to 
effect  this  quadrature. 

Vne  depth  integrals  were  evaluated  by  a  two-dimensional 
trapezoidal  rule  which  turned  out  to  be  surprisingly  accurate.  We 
?stvnate  the  calculated  values  of  ann,  to  be  within  5%  and  probably 
>  .  -  3%. 

7ne  values  of  ann,  calculated  in  this  study  are  presented  in 
Table  7.1  along  with  corresponding  values  calculated  by  Dozier  and 
Tappert.7  The  two  treatments  are  seen  to  agree  fairly  closely. 
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TABLE  7.1 


Calculation  of  a^, 
100  Hz 


n1 

a-|n,  (Penland) 

a^n, (Dozier-Tappert) 

2 

. 704E-6 

. 945E-6 

3 

. 140E-6 

. 208E-6 

4 

.525E-7 

.7601-7 

5 

. 264E-7 

.  321 E - 7 

6 

.  1 55E-7 

.  147E-7 

50  Hz 

n' 

a i n , ( Penland ) 

aln , (Dozier-Tapnert) 

2 

. 246E-6 

. 294E-6 

3 

. 548E-7 

. 710E-7 

4 

. 210E-7 

. 283E-7 

5 

. 104E-7 

. 140E-7 

6 

. 603E-8 

. 775E-8 

from 

Dozier  and  Tappert^ 

have  been  divided  by  2  to  make 

consistent  with  ours  their  definition  of  a 

nn 


Chapter  8 


Mode  coupling  induced  by  internal  waves  redistributes  the 

energy  contained  in  an  acoustic  mode  field.  Energy  confined  to  the 

SOFAR  channel  will  be  spread  over  a  much  larger  depth  at  long  ranges. 

This  scattering  causes  the  behavior  of  the  sound  field  at 

long  ranges  to  become  independent  of  source  depth.  A  limit  on  the 

depth  information  content  of  a  cw  signal  which  has  propagated  over  a 

long  distance  has  therefore  been  established. 

This  study  presents  a  theory  which  describes  this 

redistribution  of  energy.  Path  integral  descriptions  employed  by 
„  4 

Flatte  et  al.  can  be  used  to  describe  acoustic  propagation  through  a 
random  medium  at  frequencies  high  enough  to  make  normal  mode  tech¬ 
niques  invalid  or,  especially,  intractable.  However,  unlike  these 
path  integral  techniques,  this  normal  mode  description  is  valid  at  low 
frequencies  (v  £  100  Hz)  and  does  not  fall  prey  to  problems  of 
caustics  or  to  those  connected  with  describing  the  acoustic  field 
between  the  saturated  and  unsaturated  regions. 

The  scaling  techniques  introduced  by  Kohler  and 

O 

Papanicolaou  have  been  extended  to  include  azimuthal  fluctuations  in 
the  sound  speed  field.  Two  cases  are  considered.  In  one  case,  the 
azimuthal  fluctuations  in  the  acoustical  field  were  ignored  while 


retaining  ich  fluctuations  in  the  sound  speed.  Thi 
contradiction.  This  case  was  studied  to  see  if  the 
resonance  scattering  in  the  Dozier-Tappert^ ^  theory 
without  manipulating  the  spectrum  in  such  a  manner  a 
tal  isotropy,  an  average  property,  on  each  realizati 
speed.  Also,  we  wished  to  investigate  the  necessity 
assumption,  i.e.,  <AnA*,>=0  for  n*n',  in  obtaining 
Eqs.  (5.35).  In  addition,  many  concepts  are  elucida 
would  be  obscured  by  the  more  complicated  case. 

In  the  other  case,  the  azimuthal  fluctuati 
acoustic  field  are  taken  into  account.  The  results 
study  are  then  compared  with  those  of  the  previous, 
symmetric  situation. 

The  azimuthally  symmetric  case  ’ s  presente 
Chapter  5.  The  contradiction  inherent  in  a  situatic 
variations  are  permitted  in  the  sound  speed  and  not 
field  is  explicit  in  the  fact  that  the  random  pressi 
obey  an  azimuthally  symmetric  Helmholtz  equation  wtv 
somehow  "localized"  (see  the  sentence  immediately  bt 
Appendix  D). 

Contradictions  notwithstanding,  equations 

amplitude  <  Ap > ,  for  the  autocorrelation  function 

★ 

mode  amplitude  and  for  the  correlations  <  AnAn , >  b 
ferent  mode  amplitudes  (see  Eqs.  5.6,  5.10-5.13)  we 
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The  function  <An>  oscillates  with  an  exponentially 
decreasing  envelope  (Eq.  5.32).  Hence,  the  mean  acoustic  pressure 
becomes  small  at  long  ranges  more  quickly  than  would  be  implied  by 
simple  cylindrical  spreading. 

The  correlation  functions  <AnA*,>  also  oscillate  with  an 
exponentially  decreasing  envelope  (Eq.  5.34).  This  envelope  decreases 
more  rapidly  than  that  of  the  first  moment;  the  e-folding  range  is 
about  half  that  of  <An>. 

The  autocorrelation  function  < I  A  I  >  of  each  mode  is  coupled 
with  those  of  all  the  other  modes  (Eq.  5.35).  At  long  ranges  all  of 
the  functions  <i\i  >  approach  the  same  equilibrium  value;  all  depth 
information  about  the  source  is  therefore  eventually  lost. 

Because  of  the  behavior  of  <!A  |2>  and  <A  A*.>  ,  the  second 

1  n 1  n  n ' 

2 

moment  of  the  acoustic  pressure  <|p[  >  asymptotically  approaches  an 
incoherent  sum  of  contributions  from  the  the  normal  modes  (Eq.  5.39). 

The  vertical  components  of  the  average  energy  flux  <i>  decrease  with 
range  so  that  the  average  energy  flux  at  long  ranges  is  essentially 
radial  (Eq.  5.40). 

Even  without  the  assumption  of  random  phases,  the  form  of 

2 

the  coupled  equations  for  < 1  A  I  >  is  the  same  as  that  derived  by 
Dozier  and  Tappert.6  Of  course,  Kohler  and  Papanicolaou8  obtained  the 
same  form  without  assuming  random  phases  so  this  result  is  not  at  all 
surprising.  The  form  of  Eqs.  (5.32),  (5.34),  and  (5.35)  is  a 
consequence  of  approximating  the  system  by  a  Markov  diffusion  process. 

At  any  rate,  we  can  check  the  consistency  of  assuming  random  phases 


with  the  Dozier-Tappert  results  since  the  exponentially  decreasing 

envelope  of  <A„A*,>  contains  the  same  coefficients  a  as  in  the 
r  n  n'  mn 

o 

coupled  equations  for  <|A  |  >. 

Using  the  values  for  these  quantities  given  in  Table  1  of 

Ref.  7  for  a  100  Hz  signal,  we  calculate  £  (a^+a^)  =  3.8  *  10"6  m-1, 

corresponding  to  an  e-folding  length  of  265.5  km,  or  about  25 

correlation  lengths.  Of  course,  they  only  tabulated  the  results  for 

ten  modes  whereas  there  are  generally  hundreds  of  modes  in  the  deep 

(4000  m)  ocean.  Nevertheless,  the  scattering  coefficients  between 

modes  1  and  m  decrease  with  the  difference  in  mode  number  somewhat 

more  rapidly  than  |l-m|  .  Thus,  if  we  assume  a  very  large  number  of 

modes,  we  can  overestimate  £(ami+am2)  =  4.3  x  1°~6  m”\  corresponding 

m 

to  an  e-folding  length  of  231.2  km,  or  about  21  correlation  lengths. 

On  the  other  hand,  the  gravest  acoustic  mode  has  only  one  nearest 
neighbor,  so  we  would  expect  it  to  be  highly  correlated  with  mode  2. 

A  similar  calculation  of  Z^(am5+amg)  at  100  Hz  gives  an  e-folding 
length  of  about  67.2  km,  or  about  6.5  correlation  lengths.  Hence,  we 
expect  the  random  phase  approximation  in  the  Dozier-Tappert  theory  to 
be  most  valid  for  mode  numbers  of  middle  range.  Lower  acoustic 

o 

frequencies  will  have  longer  e-folding  lengths,  varying  roughly  as  w  . 

The  calculation  of  the  e-folding  lengths  was  performed  with 
scattering  coefficients  calculated  by  Dozier  and  Tappert^  using  the 
1975  Garrett-Munk  spectrum.  Their  spectrum  differed  from  the  one 
used  in  this  study  in  that  their  horizontal  wave  number  spectrum  fell 
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off  more  rapidly  with  a  than  the  more  recently  derived  spectrum  used 

in  the  analysis  presented  here.  The  1975  Garrett-Munk  spectrum  was 

32  34  35 

revised  with  subsequent  experiments  ’  '  and  we  have  used  the 

revised  version.  On  the  other  hand,  their  bandwidth  parameter  t  is 
about  twice  as  large  as  that  indicated  by  these  experiments  and  we  did 
not  obtain  scattering  coefficients  drastically  different  from  their 
values  (Table  7.1). 

From  the  form  of  the  scattering  coefficients  (Eqs.  5.30  and 
5.31)  it  is  clear  that  they  represent  contributions  from  the  internal 
wave  field  having  a  component  of  horizontal  wave  number  in  the  radial 
direction  equal  to  the  difference  in  acoustic  horizontal  wave  numbers. 
Dozier  and  Tappert  interpreted  their  scattering  coefficients  as 
representing  contributions  from  the  internal  wave  field  with  hori¬ 
zontal  wave  numbers  equal  to  the  difference  of  the  acoustic  horizontal 
wave  numbers.  However,  since  they  integrated  their  two-dimensional 
spectrum  over  transverse  wave  number,  we  may  loosely  identify  the  two 
treatments  as  having  consistent  interpretations  in  this  respect. 

At  any  rate,  one  must  always  remember  that  this  "resonance 
coupling"  is  derived  from  the  fact  that  both  the  horizontal  wave 
number  spectrum  used  here  and  the  modified  spectrum  used  by  Dozier  and 
Tappert  possess  poles  of  no  higher  order  than  one  on  the  imaginary 

axis.  These  poles  did  not  contribute  to  the  coupling  coefficients  amn 

★ 

although  the  wavelength  of  the  oscillation  of  <AnAn,>  was  indeed 
affected  by  their  presence  in  the  spectrum  (Eq.  5.31).  A  second-order 


pole  in  the  horizontal  wave  number  spectrum  would  contribute 
additional  terms  to  amn  and  a  third-order  or  higher-order  pole  could 
call  into  question  the  approximation  of  the  system  (An)  by  a  Markov 
diffusion  process. 

Because  the  system  {An}  would  be  identically  constant  in  a 
noise-free  environment  and  since  the  sound  speed  fluctuations  are 
small,  continuous,  and  have  a  reasonable  correlation  length,  it  is 
reasonable  to  expect  a  Markov  diffusion  process  to  approximate  A  with 
proper  scaling.  The  assumption  that  such  a  process  exists  itself 
places  some  restrictions  on  the  form  of  the  horizontal  wave  number 
spectrum,  but  there  is  still  a  good  deal  of  freedom  in  choosing  an 
expression  to  fit  the  measured  data.  Hence,  although  the  expressions 
for  am,,  and  bm„  can  be  trusted  to  give  accurate  numerical  results,  we 
should  remember  that  the  "resonance  coupling"  might  have  been  modified 
had  we  used  another  empirically  derived  spectrum  which  also  fits  well 
the  measured  internal  wave  data  but  has  different  analytic  properties. 

With  a  thorough  understanding  of  the  azimuthal ly  symmetric 
case,  we  now  turn  to  the  more  realistic  situation  where  azimuthal 
fluctuations  in  the  acoustic  field  are  taken  into  account.  Although 
it  is  no  longer  necessary  to  "localize"  the  Helmholtz  equation,  we  pay 
for  this  consistency  with  a  large  increase  in  complexity.  The  details 
of  this  calculation  are  presented  in  Chapter  6. 

Many  of  the  results  found  in  the  azimuthal ly  symmetric  case 
have  direct  analogies  in  the  more  complicated  study.  Each  mode 
amplitude  is  a  function  of  azimuth  for  each  realization;  each  is 


therefore  expressed  in  terms  of  a  set  of  azimuthal  modes.  The  average 
pressure  <p>  due  to  a  point  source  is  precisely  the  same  as  in  the 
azimuthally  symmetric  case.  However,  we  have  more  information:  the 
calculation  of  <p>  does  not  depend  on  an  azimuthally  symmetric 
source.  We  find  that  <p>  preserves  the  directionality  of  the  acous¬ 
tic  source  and  that  each  mode  amplitude  oscillates  in  range  with  an 
exponentially  decreasing  envelope. 

When  the  source  is  an  azimuthally  symmetric  point  source  it 

is  found  that  the  acoustic  pressure  is  a  stationary  random  process  in 

the  horizontal  plane.  The  autocorrelation  function  for  the  depth 

modes  have  exactly  the  same  behavior  as  in  the  azimuthally  symmetric 

case  provided  the  range  functions  £  <|Anq|2>  are  interpreted  as  a 

sum  over  all  the  azimuthal  modes;  energy  is  eventually  equipartitioned 

among  all  of  the  depth  modes  (compare  Eqs.  5.35  and  6.38).  As  far  as 

the  correlations  between  different  depth  modes  are  concerned,  the 

quantity  £<  AnqA*,q>  behaves  as  does  <  ApA*, >  in  the  azimuthally 

symmetric  case;  it  oscillates  in  range  with  an  exponentially 

decreasing  envelope  (Eq.  6.39).  However,  the  interference  terms  in 

<  |  P |  >  and  <  j_>  oscillate  in  q  ;  the  quantity  1^  <  AnqA*  , q  >|  must 

be  regarded  as  an  upper  bound  to  these  interference  terms.  Hence, 

azimuthal  scattering  of  energy  tends  to  slightly  enhance  phase 

randomization  so  that  the  interference  terms  are  encouraged  to  decay. 

? 

Therefore,  the  asymptotic  behavior  of  <p>,  <|p|  >,  and 
<j>  is  the  same  as  in  the  case  where  azimuthal  fluctuations  in  the 
acoustic  field  are  not  considered.  However,  this  asymptotic  behavior 


may  be  obtained  more  quickly  than  would  be  predicted  by  the  simpler 
calculation,  though  this  difference  in  range  is  probably  not 
significant.  Criteria  for  determining  the  significance  of  the 
randomization  enhancement  are  given  in  Chapter  6. 

The  result  that  energy  is  eventually  equipartitioned  among 
all  the  depth  modes  is  dependent  upon  the  neglect  of  the  continuous 
spectrum.  Unfortunately,  the  reasons  usually  cited  as  justification 
for  this  neglect  are  precisely  those  mandating  consideration  of  energy 
coupling  into  this  radiation  spectrum:  this  continuum  region  is 
attenuated  far  more  rapidly  than  are  the  discrete  depth  modes. 

O 

Analytical  work  by  Kohler  and  Papanicolaou  indicates  that 
equipartition  will  not  be  attained  if  the  continuum  is  included. 
Clearly,  further  work  is  needed  to  establish  the  rate  of  such  coupling 
via  the  internal  wave  field. 

Since  the  measure  of  volume  absorption  in  seawater  involves 

2  3 

propagation  of  an  acoustic  signal  over  very  long  ranges  (10-10  km), 

internal  waves  will  scatter  energy  out  of  the  paths  connecting  a  point 

source  with  a  point  receiver.  This  may  result  in  an  anomalously  high 

volume  "absorption."  The  theory  presented  here  should  be  modified  to 

include  the  effects  of  absorption  and  a  numerical  calculation  which 

includes  absorption  should  be  performed. 

Finally,  one  could  apply  the  methods  of  Chapter  6  to  modify 
51 

Kohler's  results  for  pulse  propagation  through  a  random  environment 


for  the  special  case  that  the  random  sound  speed  inhomogeneities  be 
due  to  internal  waves.  The  only  differences  one  would  expect  between 
the  Kohler  theory  as  it  stands  and  the  modified  Kohler  theory  would  be 
in  the  interference  terms,  which  would  be  small  at  very  long  ranges, 
and  in  the  interpretation  of  the  depth  mode  autocorrelation  functions, 
each  of  which  would  now  be  a  sum  of  azimuthal  modes. 
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Munk's  derivation  of  the  mean  sound  speed  profile  c(z)  is 

presented  here. 

Neglecting  second-order  terms,  the  fractional  sound  SDeed  gradient 
can  be  written 


l.A.l) 


Here,  T  is  temperature,  S  is  salinity,  o  is  pressure,  and  are 

constants.  Similarly,  the  fractional  density  gradient  is 


(A. 2) 


where  a  is  the  coefficient  of  thermal  expansion  and  b  is  the  coefficient 
of  saline  contraction.  The  subscript  P  denotes  "potential  gradient", 
as  explained  in  Chapter  2.  Only  the  potential  density  gradient 
contributes  to  the  stability  of  the  water  column;  only  the  potential 
sound  speed  gradient  contributes  to  sound  fluctuations  associated  with 
internal  waves.  It  is  therefore  convenient  to  introduce  the  quantity  y^, 
the  fractional  sound  speed  gradient  in  an  isentropic,  isohaline  ocean. 


Values  for  coefficients  typical  of  conditions  in  the  ocean  are  given  in 
Table  A.l.  However,  it  should  be  noted  that  a  may  vary  by  as  much  as 


TABLE  A. 1 

Coefficients  used  in  Derivation  of  Munk  Profile 


50%  between  the  sea  surface  and  the  ocean  floor.  The  value  of  a  given 
at  the  bottom  of  Table  A.l  represents  a  channel  value  for  a. 


The  relative  contributions  of  salt  and  (potential)  temperature  to 
the  stability  of  the  ocean  is  introduced  through  the  "Turner  number"  r. 


(A. 4) 


We  are  now  in  a  position  to  combine  Eqs.  (A.1-A.4)  to  obtain  an 

_  _  1  'NTT* 

expression  for  c  .  Using  the  Brunt-Vaisala  frequency. 


N2(z) 


(A. 5) 


we  find 


where 


(A. 6) 


(A. 7) 


The  sound  speed  channel  axis  is  located  where  x—  vanishes,  that  is, 

d  Z 

where 


N(z0)  ■ 


(A. 8) 


For  exponentially  stratified  oceans, 


N(z)  =  N  e"z/B 


(A. 9) 


=  i  B  In 


(A. 10) 


For  i  'surface  extrapolated"  value  of  N(0)  =  N  =  3  h  ,  B  =  1.3  km 
-2  -1 

and  Y/i  =  1.14  x  10  km”  (1  km  of  seawater  exerts  99  bars  of  nressure), 


zq  =  1.16  km 


+  as  ’ 

♦ 1 s  "^4- 


(S.ll) 


which  is  generally  about  the  same  size  as  B. 

If  we  define  the  dimensionless  variables, 


(z-zQ) 


7ByA 


the  Munk  sound  speed  profile  is  found  to  be 


c(z)  =  c(zQ)  exp[c(A+e”Li-l )] 


or,  for  small  e, 


(A. 12) 


(A. 13) 


(A. 14) 


c(z)  =  c(z0)[l  +  c(A+e”"J-l)  +  0(c2'4)] 


(A. 15) 
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We  write  our  three-dimensional  particle  velocity  due  to  internal 

4  ?5 

waves  u^  =  (u,v,w).  In  keeping  with  standard  notation,  u  is  not  the 
magnitude  of  u_. ,  but  rather  its  x  component.  Thus 


„  ,  3D. 

£}L  +  (fxu.)  +  JL  _I  =  0 

3t  1 -1  y  3x 

A  0 


1  9pi 


§?  +  (fxu.)  +  ^  TT-  pi  , 

at  -  -t  z  o0  az  i  o0 


o 


3°.  3on 


Vlj  =  o 


From  (B.2)  and  (B.3), 


3  l  32w 


3t \ay2 


~  ) 


dZ  c-.v  ' — 1 


1  9  pi  h 


°o  ?y2  \°o  3z 


o  Y-y 


(B.l) 


(B.2) 


(B.3) 


(B.4) 


(B.  5) 


(B.6) 


Similarly 


_L  ( .  _L  ly.  ]  +  (fxu. ) 
3tU2  3z3x/  3X^  1  ‘ 


-1-  -L  (fxu.) 
3z  3x  -  -1\ 


+  — - 


32Di  ,  1  30^ 
— ?  ^  3z 


o  3x 


_  5. 

o„ 


(B.7) 


Now 


1  "Pi  _  3u 


PQ  3x  3t 


(B.8) 


and 


(ixy-i) 

pQ  ay  3t  y 


(3.9) 


Using  the  incompressibility  condition  (B.5), 


—  (72w)  +  72(fxui- )  -  [7  •  (fxUj)] 


(A*  3? 


°o  xlx7  iy7/  O  3z 


[_3_  3w  _ 

3t  3z 


(fxu. ) 


(3.10) 


We  now  define  the  vorticity  l. 


£  =  7*u. 


(B.ll) 
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V 

‘  (1*^)  = 

-f  • 

t 

and 

V2(fxU  ) 

z 

=  - 

&<!•£> 

+  V 

•  £z*(v 

This  last 

term  involves 

zx[yx(fxu 

,)]  = 

:  x(f  •  V)v 

-  y(f  •  v)u 

Thus 

&  (?2*) 

+  (\W  •  fy  4) 

[fxfe  + 

Z7i 

lv 

-y  / 


■■o  -•  L  •  '■ 


2  2  2 
For  N  >>qV<;  , 


K 


Thus,  differentiating  (B.15)  with  respect  to  time  and  employing  :he 
incompressibi 1 i ty  condition  (B.5), 


4(v2«)  *n2(Z)  (£*£)»  +i- 

3t“  \3x“  3y  / 


N2(z) 


JL_  M  +  f  .  + 

at2  3z  1  9t 


3t3Z 


(fxu^). 


The  vorticity  equation,  which  is  the  curl  of  equations  (B.l), 
(B.2),  and  (B.3)  added  together,  is 


3 1 

at 


V  x 


(fxu. )  +  7 


( 


and,  since 


we  have 


V  x  (fxu^ )  =  -  (f  •  V)Ujj 


(I 


->y  (V2w)  +  N2(z)  +  -^Aw  +  (f'V)2W 

at  \3x  3y  / 


d  3w  ,  -  / 

n  Pi9  p\ 

Li?w 

V  x  -  z 

Do  / 

+  f  •  (f -V)^  +  4^7  (fxuj_  .  ( 


APPENDIX  C 


Much  of  the  following  discussion  of  Markov  processes  comes  from 

q 

Ludwig  Arnold's  marvelously  lucid  book  and  from  lecture  notes  given  by 
57 

Werner  Horsthemke.  More  mathematical  detail  than  is  given  in  the 
discussion  here  may  be  found  in  these  references. 

The  Markov  property  of  stochastic  dynamic  systems  says  that  "if 
the  state  of  a  system  at  a  particular  time  s  (the  present)  is  known, 
additional  information  regarding  the  behavior  of  the  system  at  times  t<s 
(the  past)  has  no  effect  on  our  knowledge  of  the  orobable  development 
of  the  system  at  t>s  (the  future)". 

As  an  illustration,  let  us  consider  a  discrete  system,  a  machine 
which  at  the  touch  of  a  button  chooses  at  random  an  integer  from  the 
closed  interval  [0,9].  This  machine's  inner  workings  will  act  so  that 
the  Drobability  of  getting  any  particular  integer  at  the  next  push  of 
the  button  depends  only  on  what  the  last  number  was,  and  not  on  any  of 
the  numbers  given  before  that.  For  example,  say  I  have  pushed  the 
button  three  times  at  times  t^ ,  t9,  and  t^  with  the  resulting  sequence 
(9,4,6).  Now,  knowing  that  the  last  value  was  6,  the  probability  of 
getting  a  7  is,  say,  30%.  That  is, 

P(7,t4l6,t3)  =  0.3  .  (C.l) 

Let  us  also  say  that  (C.l)  would  have  also  been  true  if  my  previous 
trials  had  given  (1,2,6),  (3,5,6)  or  any  sequence  ending  in  6.  Then, 
our  machine  has  the  Markov  property. 


As  an  example  of  a  non-Markovian  system,  let  us  say  that  in  order 
to  get  a  7  after  a  6,  it  is  necessary  to  qet  a  5  before  the  6  and, 
having  gotten  this  5-6  combination,  I  still  have  a  30%  chance  of  getting 
c  7.  Then 

P(7,t4l(6,t3),(5,t2))  -  0.3  (C.2) 

and,  for  example, 

P(7,t4|(6,t3),(8,t2))  =  0.0  .  (C. 3) 


It  is  pretty  hard  to  get  a  7  after  a  6  in  this  system!  But  the  lack  of 
the  Markov  property  is  clear;  knowledge  of  the  past  improves  our 
prediction  for  the  future.  (Note;  We  do  not  consider  the  so-called 
"higher-order  Markovian  systems"  to  be  true  Markov  processes.) 

Notice  that  it  is  not  necessary  for  each  trial  to  be  totally 
(statistically)  independent  of  all  of  the  others  for  the  system  to  be 
Markov.  The  Markov  property  simply  states  that  given  a  set  of  trials 
and  outcomes,  ((xn,tn) , (xn_1  ,tf)_1 ) , . . .  ,(x]  ,t] ))  , 


P(xn+1  ’Vl  i  ^xnftn^  ’^xn-l  ,t:n-l )  *  •  •  •  ’<xi  ))  =  p(xn+l » tn+l  I  xn ’  tn^ 


(C.4) 


P(xn+.|  ,tn+-|  |xn,tn)  is  called  the  "transition  probability". 

The  Markov  property  can  also  be  true  for  continuous  systems. 

Since  we  are  talking  about  stochastic  systems,  we  must  be  more  specific 
about  what  is  meant  by  "continuous",  and  that  necessitates  the  concept 
of  a  "transition  probability  density".  Consider  a  one-dimensional 
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system.  If  we  know  the  transition  probability  P(B,t|x,s)  that  the 
system  will  be  in  some  interval  B  at  time  t  if  it  is  in  state  x  at 
time  s<t,  the  "transition  probability  density"  p(b,t|x,s)(if  it  exists) 
is  defined  by 


P(B 


,t|x,s)  =  / 
B 


db  p(b.r ! x  ,s) 


(C.5) 


We  shall  be  concerned  with  systems  that  are  continuous  "almost  surely": 


lim 

t+s 


dy  p(y,t|x,s)  =  0 


(C.6) 


y-x  |>£ 


That  is,  as  the  length  of  time  between  measurements  decreases  to  zero, 
the  probability  of  finding  the  system  in  states  separated  by  some 
difference  e  goes  to  zero  more  quickly  than  (t-s),  no  matter  how  small  e 
is.  We  are  finally  in  a  position  to  define  a  Markov  diffusion  process. 

A  Markov  diffusion  process  is  an  almost  surely  continuous  Markov  process 
where  the  "drift"  f(x,s)  and  "diffusion"  g  (x,s),  given  by 


\--\i 

v-'-j 


-I!. 


exist.  In  many  cases  of  interest,  such  as  the  cases  we  consider,  the 
cutoff  e  can  be  extended  to  the  entire  range  of  y. 


The  importance  of  these  processes  lies  in  that,  since  our  system 
is  a  Markov  diffusion  process,  we  can  immediately  write  a  (backward) 
propagation  equation  for  the  moments.  Often,  as  in  the  cases  we 
shall  consider,  we  can  immediately  write  a  (forward)  propagation 
equation  for  the  transition  probability  density,  or  "Fokker-Planck" 
equation: 

dp'^dt^X~,'~  =  '  37  n(.v,t|x,s)  +  \~2  Q2( V »t)  p(y,t|x,s) 

9  V 

(C.9) 

Now  the  reason  our  system  is  called  a  " diffusion  process"  is  clear  from 
the  form  of  the  Fokker-Planck  equation,  which  is  that  characteristic  of 
deterministic  diffusion  equations. 

More  general  treatment  of  these  concepts  can  be  found  in  excellent 
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books  by  Stratonovich,  Arnold,  or  Horsthemke  and  Lefever.  In 

particular,  we  shall  need  the  Fokker-Planck  equation  analogous  to 

Eq.  (C.9)  for  a  system  involving  several  coupled  processes  (the  acoustic 

g 

normal  modes).  Kohler  and  Papanicolaou  have  shown  the  usefulness  of 
the  theory  in  acoustics  and  we  shall  apply  it  to  the  special  case  of 
acoustic  propagation  through  a  field  of  internal  wave-induced  sound 
speed  fluctuations. 


APPENDIX  D 


In  terms  of  the  components  of  A,  Eq.  (5.28)  becomes 


dp(A,A*;a|A0,A^;a0) 


do 


ii  p*H 


1  im 

Tfoo 


nm  'a 


n'm' 


fsTT  Vn'm  ‘  MI  ITT  AnAn'l 

L  m  m  m  J 


mm  - 

+  <Vnm(t)VnV(s)>  AnVm  *  3^  jff-  AX'. 


+  < V  (t)V*,  ,  (s) >  -4  r  A  A* 
nmv  1  n'm''*1  37T  afl*7 
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<V™(t>Vm'<s»  5**3)rr*nV}  | 

mm  J 


;oQ)  ,  ( 


with  p(AfA*;a0|A0,Aj;o0)  =  6(A-A0)5{A*-Aj) 


To  get  a  propagation  equation  for  <A^>  ,  we  multiply  (D.1)  by  A| 


and  integrate  over  A  and  A*.  Integration  by  parts  gives  us 


d<Ak> 


do 


.  ro  +T  r  S 

S-Z11wf  °  ds  dt<V_(t)V„lml(s)>  <A„>5„lm5,m, 


nm '  n  m 


n 


From  Eq.  ( b - 1 4 ) 


<V„m(t)V  .  ,(s)>  = 
rrnr  '  n  m 


rH  th 

i  -I, 


k  k  ,k,  o 


n  m  n  m 


P  ( z )  r  ( z 1 ) 


<60(1,2)60(5,2'  )> 
c3(2)c3(2'  ) 


i'(k  -k  )t  l ( k  1  - k  ,)s  .  (0.3) 

^mn'm  n 
e  e 


If  s  and  t  are  in  the  same  direction,  i.e.,  the  line  connecting  source 
and  receiver, 


<6c(t,2).6c(s,2_'J>  =  8*  <^2(o)>  N.2(2|N^(2^ 


c  ( 2 )  c  ( 2 1  ) 


f  tt/2  rN(n)  r 00  a 

d0  dQ  d,  -y— o 

J0  Jf  J- 00  a  +  a* 


J/2 


cos(n^x(a))  eiacose(t-s) 


(D.4) 


0/2 


/  2  2 

where  n  =  ^  (z+2'),  ^  =  [2-2'|,  X  =  (  ~yn— — g-— )  and  the  other 

quantities  are  defined  in  Chapter  3. 


Let  Ak,=k  -k  ,  Ak0=k  ,-k  , .  Consider 
I  m  n  cm  n 


ro  +t  fs  i (acose+Ak, )t  i (Ak?-acos0)s 


n 


1(Ak,*Ak2)00  (ei(Ak,*Ak2)T  J  1 


/  i  (Ak2-occos0)T  \ 
i  ( Ak]  +Ak2  )oQ  (acosQ+Ak, )  (Ak0-a.cos6; 


(0.5) 


Now  we  consider  the  a  integration.  Because  of  the  operation  lim  1/T, 

Tfoo 

only  terms  which  vary  at  least  linearly  with  T  will  survive.  Hence, 
we  exhibit  only  results  of  the  a  integration  which  meet  that  criterion. 

From  Eq.  (D.5),  I  consists  of  two  terms.  Consider  the  first 
a  integration. 


(^V4k2)T  .  ,)  i(Ak,*Ak2)a 


cos  (a 


[Ak1+Ak2; 


+a^  (acos9+Ak.|) 


(D.6) 


The  only  way  that  I  can  be  prooortional  to  T  is  if  Ak^-Ak^.  Thus 


[1  =  -  iT( 
•/-a 


cos [a  ?  x(n) ] 


^  (•  * 


«5(Ak1+Ak2)  •  (D.7) 


Note  that  the  pole  at  a  =  -  ■-  appears  to  lie  on  the  real  axis. 

cost) 

However,  the  quantity  Ak-j  will  have  a  small  imaginary  part  if  the  ocean 
is  slightly  absorptive.  We  will  use  this  imaginary  part  to  lift  the 
pole  off  the  real  axis  and  then,  after  evaluating  the  integral,  let 
this  imaginary  part  go  to  zero.  Since  Ak^k^k^  Ak-,  will  lie  in  the 


fourth  quadrant  of  the  complex  plane  for  n> m  and  in  the  second  quadrant 

A  A 

of  the  complex  plane  for  n<m.  Thus, 


I-|  =  -6(Ak.|+Ak2)TrT 


-ot  .ft  i  Ak,  X  4/cose 

iAk-|  e  +  °t  e  cose  .  . 

- - - - - s -  ,  n<m  ,  (0.3) 

/.  _ 


(c£cos  9+AkJ) 


'(Ak,  ^  /’-V  ?'“»« 


I-j  =  -6(Ak-|+Ak2)^T 


cose  l  ,  n>m  . 


(a  cos^G+Akf) 

★  I 


The  second  a  integration  is 


h  *  'e 


i(Ak1+Ak2)o0  r<=°  dQ  C0S(Q^X)(el(Ak2  aCOS0)T  _]) 

-L  cos2e  („  2A  2\(  .  Akl  V  Ak2  \ 
r  +vr  cose/r  '  cose  j 


.  (0.10) 


The  only  way  (D.10)  can  be  proportional  to  T  is  for  the  inteqral  to 
contain  a  double  pole;  i.e.,  Ak^=-Ak2< 

flow  we  see  a  problem.  In  Eq.  (0.10)  the  quantity  X(O)  varies  from 

zero  to  infinity  before  we  take  our  limit  T-*».  We  would  like  to 

evaluate  I2  by  contour  integration  and  the  relative  sizes  of  Tcose 

and  ^X  will  make  a  difference  in  the  way  we  close  our  contours.  Thus, 

we  must  consider  Eq.  (0.4)  in  order  to  estimate  the  Q  integration 

between  \  cose  and  X  -*on  # 
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Let  us  define  a  small  quantity  5  such  that  Tcos-  •  fcr 
O<0<tt/2-6  and  ^X>Tcos8  for  tt/2-6<9<tt/2. 

Then,  dropping  the  terms  independent  of  T, 
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V  cos*1  ' 


:cos  •'*+ 
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(0.17) 


1-Je  may  therefore  neglect  the  remainder  term  for  very  small 
Physically,  we  expect  this.  The  integration  procedure  combined  with 
the  limit  on  T  picks  out  of  the  internal  wave  soectrum  those  radial 
wave  number  components  (components  in  the  direction  of  acoustic 
propagation)  equal  to  (k^-k^l.  By  neglecting  the  term  we  say 

that  internal  waves  traveling  in  a  direction  orthogonal  to  the  direction 


of  acoustic  propagation  do  not  cause  scattering  of  energy  among  the 


acoustic  depth  modes.  That  is,  the  probability  of  finding  internal 
waves  with  infinite  wave  number  is  zero. 

If  we  impose  the  condition  that  acoustic  eigenvalues  k  along  with 
their  sums  and  differences  exist  uniquely,  then 


^(Aki+Ak2)  -  6nm.6mn, 
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Equation  (D.19)  has  the  solution 


<Ak(^)>  =  exp 
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where 

\  ^amk-1^mk^ 

m 


Note  that  even  with  the  oscillatory  terms  amk>0  and  is  symmetric 
and  k  while  b k  is  antisymmetric  in  m  and  k.  Note  also  that  Ak 
independent  of  c  . 


is 


(D.21) 
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APPENDIX  E 

• 

Consider 

1 

• 

2  sin[J  (Aq+Aq')] 
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Since  Aq,  Aq 1  are 

integers,  the  expression  (E.l)  is  equal  to 

• 

^  3  sin[^  (Aq+Aq' )] 
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.  • 

9  In  the  second  integral,  let  <p-*--<p,  so 

_  • 
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If  we  perform  similar  manipulations  with  the  q>'  inteqral. 
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(E. 

Since  P(a,Q;n)  cos(a^X)is  even  in  a,  the  integration  over  a  will  have 
nonzero  contribution  only  if  Aq  and  Aq'  are  both  even  or  both  odd. 
However,  if  they  are  not  also  equal  in  magnitude  and  ODposite  in  sian, 
the  sine  term  will  vanish.  Thus. 
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APPENDIX  F 


Consider  the  expression 
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The  j=l  term  is 
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as  long  as  A and  AQ2  are  ' 
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Again  we  see  the  operator  evaluated  in  Appendix  D  and 
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T  Z,  J-i  nq  k  j  p  ,q-p+q  m  k  p  j  mn  f 


1  mpnq  1 
m'p'n'q1 


"tS  (amk"ibmk^  <AkjAk‘ j 


The  second  term  is  the  complex  conjugate  of  the  first  with  an  exchange 
of  indices  ( k,j )+-►( k1  ,j' ): 


T  ^<AnqAkj>{mn'{pql6m'k'Vj'V.q-P+ql 


m'p'n'q* 


"  tX  (amk'+  ^mk'5  <AkjAk'j,> 


Calculation  of  the  third  term  is  a  bit  trickier: 


mpnq  < AnqAn'q (>  ^mk^pj^m' k ' V j ' 6p ’ ,q ' -q+p 


m'p'n'q' 


I-"*  •**  *  *•’ 


r-*  h  2  7  fCT  cc"  rM  rH 

Z  ~TT  (a)  V<c(0)>  d*'  dz  dz'e(2)e(; 

™  e  tt  Jc  Jo  Jo  Jo 


<Pn(z)<pk(z)<pn,  (z*  )<Pk,  (z*  )N2(z)N2(z‘  )  rN(n)  /*  » 


v/knkkkn,kk'  c^2)^2' )N(n) 


I  dQ  |  da  cos  (a  ^X) 


(P.2-f2)1/2  ?. 


fVI  ?.  f  VT  ds  I" 5  dt  ,<  ( J-q  V  .-1  ( i  ' ■  -q  • '  * 

sr  a  +  a*  Ja  Ja 


"^kk“kn^s  _1 (kk'"kn' ^  ia(scos^-tcos^' ) 


<  AnqAn'q' >  6 j ' - j .q'-q 


A  similar  term  is  obtained  for  the  fourth  term.  It  is  here  that 
excitation  of  the  azimuthal  modes  is  explicit.  Even  if  only  the 
azimuthally  symmetric  mode  is  initially  nonzero,  all  of  the  other 
azimuthal  modes  are  excited  as  long  as  j-j1  is  zero.  This  null 
difference  between  j  and  j'  is  maintained  throughout  the  entire  field 

p 

of  propagation  and  the  expression  for  <!d|  >,  for  example,  becomes 
.2,.  2  2^  <"AkiAk'i'> 

<|p|  >=u  p  J  — - *J-  k  J  -  dk(z))k,(z') 


k*  vVk' 


|l  F(  k.  -k.  ,  )r  +  dfX  -  1  fll 


That  is,  for  an  azimuthally  symmetric  source,  the  different  azimuthal 
modes  are  uncorrelated  to  first  order  in  e. 

After  much  manipulation,  we  find  that 


d<Vk'j> 


"  '  hi"*'1  +  1  ( bmk •  )]<AkjAk,j> 


KV',V)  •  (r-7) 


where  the  sums  over  m  and  m’  are  restricted  such  that  I  kk-kml  ^  I  kk'  ~km' 
and  that  (kk>-km«)  has  the  same  siqn  as  (kk~km).  The  quantities 

dSkn,  •  k  •  are  <|1',en  ^ 


^  ...  =  /U\2  M  f  sr2[ 


H  r  H 


mkm'k' 


(|f  V<^°»  _[  d2£ 


1(z)d>m.(z,)4>k(z)d>|<.  (z* ) 


kmkm,kkkk' 


fN(n)  fo2.f2l1/2 
p(z)o(z* )  /  dsii— V - 


■  .N2.Lz)N2[z') 
r(z)c2(z')N(n) 


x  I  d<A:os<£cos[(  j-q)^]  cos  (j-q)cos 


x  cos 


kk'"km' 


V-ip-  ?*! 

(a2cosV+  (kk.-km,)2) 


.  (G.8) 


.  •  -  •  .v. 


The  expression  (G.7)  is  prohibitively  difficult  to  calculate  since 
the  terms  would  have  to  be  evaluated  M4(2L+1)2  times  -n  order 

to  specify  the  propagation  of  each  term!  However,  much 

information  can  be  obtained  about  energy  transfer  without  actually 
solving  for  each  <AkjA*.  ■>  .  From  Eq.  (G.6),  we  see  that  < | p | 2 > 
(and  also  <j_>)  contains  terms 


"k  ■ 1  <V2> 


If  we  sum  Eq.  (G.7)  over  j  we  find  that 


uWj^  p 

"d5~  “  amk(Vwk} 


(G. 10) 


if  k=k '  and 


lr£  <AkA'j> 

j 


7  S  [(amktamk'>  +  ,'(bmk'-bmk>]^<AkjAk'j> 


for  k^k ' 
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